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ABSTRACT 

The  airborne  transmission  loss  across  resilient  mounting  systems  is  investigated 
by  developing  an  analytic  model  of  the  transmission  of  vibration  from  a  point-driven, 
simply  supported  plate  to  a  parallel  infinite  plate  through  the  air  separating  the  plates. 
Predictions  of  airborne  transmission  losses  are  compared  with  transmission  losses 
through  resilient  mounts  to  quantify  the  role  of  the  airborne  path  in  transmission  losses 
through  well  designed,  single-stage  resilient  mounting  systems.  The  sensitivity  of  the 
airborne  transmission  losses  to  separation  distance,  thickness  of  the  receiving  plate,  size 
of  the  driven  plate  and  damping  in  the  plates  is  presented. 
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NOMENCLATURE 


c  speed  of  acoustic  wave  propagation 

D  bending  rigidity 

E  modulus  of  elasticity 

f  frequency  (Hz) 

F  input  force 

G(r/rs)  Green's  function 

h  thickness  of  the  plate 
i  square  root  of -1 

k  acoustic  wavenumber  or  stiffness 

kp  free  bending  wavenumber  in  the  infinite  plate 

Lg  length  of  plate  A  in  the  x  direction 

Ly  length  of  plate  A  in  the  y  direction 

xQ,  yQ  drive  point  on  plate  A 

z  distance  between  plates 
w  plate  response 
Y  ratio  of  specific  heats  =  1.4  for  air 
T|  damping  factor 

p  viscosity  of  fluid 

v  Poisson's  ratio 

p  density 

to  frequency  (rad/s) 

°>o  resonance  frequency 
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Resilient  mounts  are  often  used  to  reduce  the  propagation  of  unwanted  vibration 
from  a  machine  to  its  foundation  or  to  reduce  the  transmission  of  motion  of  a  foundation 
to  vibration  sensitive  equipment.  Most  research  on  and  development  of  resilient 
mounting  systems  has  focused  on  the  structure  borne  paths  through  the  mounts.  The 
often  neglected  airborne  path  between  the  vibrating  machine  and  foundation  is 
considered  in  this  thesis  by  developing  an  analytic  model  for  the  airborne  propagation  of 
vibrations  from  a  finite  elastic  plate  to  a  parallel  infinite  plate.  In  the  model,  the  unsteady 
pressures  generated  by  vibration  of  the  upper  plate  are  propagated  to  the  surface  of  the 
infinite  plate  where  the  incident  unsteady  pressures  excite  the  infinite  plate  into  vibration. 
The  airborne  transmission  loss  is  then  computed  as  the  ratio  of  the  amplitudes  of  the 
vibrations  of  the  two  plates. 

A  simple  diagram  of  a  mounting  system  using  four  rubber  mounts  is  shown  in 
Figure  1.1.  Using  four-pole  parameters  as  developed  by  Molloy  (1),  the  transmissibility 
of  vibration  from  a  vibrating  rigid  body  foundation  with  uniform  velocity  to  an 
unconstrained  mass  is  given  by 

Vj  im  I  +  Op2 

T  -  Vj  * .  C  2  2  ’ 

»®  m  +  °>o  -  ® 

where  the  undamped  natural  frequency  is  given  by 

“o  = 
and 

C  =  4  times  the  damping  factor  for  one  mount 
k  =  4  times  the  stiffness  of  one  mount 


1.1 


1.2 


rigid  body  (lumped  mass) 


Figure  1.1  Simple  diagram  of  a  mounting  system 
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m  =  the  mass  of  the  mounted  item. 

Using  equation  1.1,  the  vibration  transmission  losses  for  propagation  through  mounts 
with  different  natural  frequencies  are  plotted  in  Figure  1.2.  According  to  equation  1.1, 
the  transmission  loss  increases  indefinitely  as  one  over  the  frequency  squared  at 
frequencies  much  greater  than  the  natural  frequency.  However,  measurements  on  single 
stage  mounting  systems  have  rarely  shown  transmission  losses  as  high  as  the  predicted 
losses  given  in  Figure  1.2. 

There  are  several  reasons  why  the  transmission  losses  predicted  by  equation  1.1 
are  not  achieved  in  reality.  Some  of  these  are: 

1.  Structurebome  flanking  transmission  paths  through  mechanical  connections  to 
the  mounted  machinery,  such  as  piping,  electrical  connections  and  duct  work. 

2.  Wave  effects  in  the  mounts  at  high  frequencies  where  resonances  inside  the 
mounts  occur  increase  the  transmissibility  through  the  mounts.  This  effect 
has  been  previously  studied;  see,  for  example,  Snowdon  (2). 

3.  Airborne  transmission  from  the  vibrating  structure  across  the  mounts  to  the 
vibration  sensitive  structure  or  equipment 

Because  the  last  effect  is  often  neglected  in  mount  design  and  because  the  airborne  path 
may  be  a  significant  factor  in  the  transmissibility  for  single  stage  mounts  at  high 
frequencies,  it  will  be  the  subject  of  this  thesis. 

Heckl  (3)  measured  the  vibration  transmission  loss  from  a  wall  to  a  plate  parallel 
to  the  wall  with  and  without  mounts.  The  transmission  loss  was  first  measured  with  the 
plate  mounted  on  rubber  mounts  having  a  fundamental  resonance  of  40  Hz.  The  mounts 
were  completely  removed  and  the  transmission  loss  measurements  repeated  with  the 
same  distance  between  the  plate  and  the  wall.  Figure  1.3  shows  the  transmission  losses 


Transmission  Loss  (dB) 


38  70  120  mm 


63  125  250  500  1000  2000  4000 


Frequency  (Hz) 


Figure  1.3  Transmission  loss  for  a  plate  mounted  on  a  wall  with  and  without  rubber 

^  Jnounts  (3) 
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between  the  plate  and  the  wall  with  and  without  the  mounts.  Above  100  Hz,  the 
vibration  transmission  losses  with  the  mounts  are  nearly  equal  to  the  losses  without  the 
mounts.  This  suggests  that  energy  transmitted  between  the  plate  and  the  wall  via  the 
airborne  path  is  nearly  equal  to  the  energy  transmitted  from  the  plate  to  the  wall  via  the 
structure  borne  path  through  the  mounts. 

Ungar  (4)  treated  the  air  beneath  the  isolated  equipment  as  adding  stiffness  to  the 
stiffness  of  the  resilient  mounts,  thus  increasing  the  stiffness  of  the  mounting  system. 
The  resonance  frequency  of  the  entire  system  including  the  resilient  mounts  and  the 
entrapped  air  is: 

f„  +  1-3 

- 

In  the  absence  of  air  the  resonance  frequency  fs  of  the  mounting  system  is: 

fs  -(£)Vk*n  14 

where 

kg  =  the  static  stiffness  of  the  mounts 

m  =  the  mass  of  the  mounted  item. 

Without  the  mounts,  the  resonance  frequency  fa  of  the  entrapped  air  is: 

fa=(^)  15 

*a  =  Wd 

where 

kg  =  the  stiffness  of  the  air 

y  =  ratio  of  specific  heats  =  1.4  for  air 
pQ  =  the  ambient  air  pressure 

d  =  the  air  gap  thickness 


Transmission  Loss  (dB) 
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Figure  1.4  shows  the  classical  transmission  loss  curve  without  the  air  (curve  A) 
and  the  transmission  loss  with  the  combined  stiffness  of  the  mounts  and  entrapped  air 
(curve  B).  These  curves  indicate  that  the  airborne  path  may  be  significant  in  the 
transmission  across  resilient  mounts. 

The  results  in  Figure  1.4  are  only  valid  when  the  air  is  entrapped.  If  the  space 
between  the  mounted  equipment  and  its  foundation  is  vented,  then  the  inertia  and 
viscosity  of  the  air  must  be  taken  into  account  when  calculating  the  transmission  loss.  If 
the  air  is  allowed  to  escape  freely  from  under  the  edges  of  the  mounted  equipment,  the 
contribution  by  the  air  to  the  total  stiffness  may  be  less  at  low  frequencies. 

The  above  analysis  does  not  apply  at  higher  frequencies  because  it  assumes  that 
the  base  of  the  isolated  equipment  moves  as  a  rigid  body.  At  higher  frequencies  the 
equipment  and  foundation  no  longer  behave  as  rigid  bodies. 

To  include  the  effects  of  the  venting  of  air  in  between  the  mounted  equipment 
and  its  foundation,  and  the  elastic  vibratory  motion  of  the  mounted  equipment  and 
foundation,  an  analytic  model  for  the  airborne  transmission  between  parallel  elastic 
plates  is  developed  in  this  thesis.  In  Chapter  2,  the  mathematical  model  for  the  response 
of  a  finite  plate  to  a  point  excitation  and  the  airborne  transmission  loss  from  the  finite 
vibrating  plate  to  the  vibration  response  of  a  parallel  infinite  elastic  plate  is  derived. 
Predictions  of  the  airborne  transmission  loss  are  presented  in  Chapter  3  and  compared  to 
predicted  transmission  losses  through  resilient  mounts.  Also  the  dependence  of  the 
airborne  transmission  losses  on  model  parameters,  such  as  the  separation  distance 
between  the  plates,  thickness  of  the  infinite  plate,  size  of  the  finite  plate,  frequency  and 
damping  in  the  plates  are  presented.  Conclusions  and  recommendations  for  additional 
research  are  given  in  Chapter  4. 


Chapter  2 

ANALYTIC  MODEL 
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2.1.  Introduction 

In  the  analytic  model  of  the  airborne  transmission  from  a  vibrating  piece  of 
machinery  to  its  supporting  structure,  derived  here,  the  foundation  (receiver  structure)  is 
modeled  by  an  infinite  flat  plate  (plate  B)  and  the  machinery  (source  structure)  is 
modeled  as  a  finite  flat  plate  (plate  A)  parallel  to  the  surface  of  the  receiving  structure  as 
shown  in  figure  2.1.  In  developing  this  model,  the  following  steps  are  taken. 

1.  The  solution  for  the  response  of  the  top  simply  supported  plate  (plate  A)  to  a 
point  excitation  is  derived  using  classical  modal  expansion. 

2.  Euler's  equation  is  used  to  relate  the  velocity  of  plate  A  to  the  unsteady 
pressure  on  the  surface  of  plate  A  facing  the  infinite  plate  (plate  B). 

3.  Using  the  wave  equation,  the  pressure  on  the  surface  of  plate  A  is  propagated 
to  the  surface  of  plate  B. 

4.  Fourier  transforms  are  used  to  transform  the  pressure  on  the  surface  of  plate 
B  into  wavenumber  space. 

5.  Using  the  Green's  function  for  plate  B  in  wavenumber  space,  the  solution  for 
the  response  of  plate  B  to  the  pressure  on  its  surface  is  obtained  as  the  product 
of  the  Green's  function  and  the  pressure  on  the  surface  of  plate  B. 

6.  The  inverse  Fourier  transform  is  taken  to  obtain  the  spatial  distribution  of  the 
response  of  plate  B. 

7.  Predictions  of  the  transmission  loss  are  derived  by  taking  the  ratio  of  the 
averages  of  the  vibration  responses  computed  for  plates  A  and  B. 


10 


Figure  11  Infinite  flat  plate  with  a  finite  flat  plate  parallel  to  it 


2.2.  Assumptions 


One  of  the  assumptions  in  the  development  of  the  analytic  model  is  that  the 
impedances  of  the  plates  are  much  larger  than  the  acoustic  impedance  of  the  air.  With 
this  assumption,  a)  the  effect  of  air  loading  on  the  vibration  of  plate  A  will  be  neglected, 
b)  the  response  of  plate  A  to  pressure  radiated  by  plate  B  back  to  plate  A  will  be 
neglected,  and  c)  the  pressure  on  the  surface  of  plate  B  will  be  the  "blocked"  pressure, 
i.e.  the  motion  of  plate  B  will  have  little  impact  on  the  pressure  on  the  surface  of  plate  B. 
Therefore,  with  the  impedances  of  the  plates  much  larger  than  the  acoustic  impedance, 
the  solution  for  plate  vibrations  decouple  from  the  air  pressures.  Also  standing  waves  in 
between  the  two  plates  will  not  be  included  in  the  model. 

2.3.  Mathematics 

2.3. 1  Vibration  Response  of  Plate  A 

For  Plate  A,  simply  supported,  and  driven  by  a  point  source  located  at  xQ  and  yD 
as  shown  in  figure  2.1,  the  plate  equation  is 

Dav4 wa  *  Paha0)2wa  =  F&(x-Xo)5(y-yo)  »  21 


where  the  bending  rigidity  of  plate  A  is  given  by 
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where  pa  is  the  density  of  plate  A, 
ha  is  the  thickness  of  plate  A, 

Ea  is  the  modulus  of  elasticity  of  plate  A,  and 
v  is  Poissons  ratio. 

For  a  simply  supported  plate  the  boundary  conditions  are: 


and 


Wa 

-  o'! 

a2wa  > 

7? 

=  0J 

W  a 

-  o'! 

a2wa  _ ) 

dy2 

=  0j 

for  x  =  ±  1^/2 


With  opposing  edges  simply  supported,  the  solution  is  separable,  so  that 
wa(x,y)  =  X(x)Y(y)  . 

To  satisfy  the  boundary  conditions  given  by  equations  2.3, 

E(2m+1)  n  x  1 

Lx  J 

f(2n+l)  Jt  y  1 

Y(y)  -  cos  L  Ly  J  . 

The  solution  can  then  be  expressed  as 


oo  oo 

wa(x,y)  -  X  X 
m=0  n=0 


Amn  cos  [ 


(2m+l)  7t  x  i  r(2n+l)  it  y 


]  cos  P 


Hr 


] 


2.3 


2.4 


2.5 


.  2.6 
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The  resonant  frequencies  co^  are  obtained  by  using  equation  2.6  in  the  homogeneous 


equation  2.1, 

?  Da  r  r(2m+l)jc-|2  r(2n+l)7r|2i  2 

“mn  +  L~J  1  ' 

Now  using  equation  2.6  in  equation  2.1,  and  operating  with 
Lx/2  Ly/2 

f  f  cos  p,*|)«y]dxdy 

-Lx/1  -Ly/2  ^  J 

and  using  orthogonality  yields 

r(2m+l)7Cxol  r(2n+l)jty0-| 

4  Fcos[  ~~  L^— J  cos[  ~  L,  J 

Am"=  Mp[“mn2-“21 

where  Mp  —  p^ha  Ly , 

Using  equation  2.8  in  equation  2.6, 


2.7 


2.8 


wa(x,y)  = 


4F 

Mp 


oo  oo 


II 

m=0  n=0 


COS 


[ 


(2m+l)Jtxr 


1  p(2n+l)jty0-j  r(2m+l)jtxi  f(2n+l)7tyi 

JcosL~- ^  JcosL  l:  JcosL"u  J 


[  (Dmn2  -  CO2  ] 


Jx 


2.9 


Taking  the  Fourier  transform  of  equation  2.9 


oo  oo  f(2m+l)nxo1  r(2n+l)7ty0i 

4FV  r  cosL  L.  J«>sL  Ly  J 

wa(kxJcy)  =  M^2-f  -L  [  o)mn2  -  co2  ]  Imn(kx»ky)  210 

rm*0  n=0  1  mn  1 


where 


00  00  r(2m+l)7cx-|  r(2n+l)jcyi  _:i,  Y  _:l.  v 

^mn^x’^y)  =  /  JC0SL  Jcos|_  Je  x  e  y  dxdy  .  2.11 


ImnOcx’ky)  is  evaluated  in  Appendix  A. 

2.3.2.  Acoustic  Field  Between  Plates 

The  unsteady  pressure  in  the  air  between  plate  A  and  plate  B  satisfies  the  wave 
equation  which,  for  harmonic  time  dependence  e  "10)t,  becomes  the  Helmholtz  equation, 

V2p  +  k2p  =  0  2.12 


00 

where  k  =  ~  is  the  acoustic  wavenumber. 

Taking  the  Fourier  transform  of  equation  2.12  yields 


2  2  2  ® 

(k  -  kx  -  ky  +  ^2)  p  (kXJky;z)  =  0 


2.13 


The  solution  to  equation  2.13  is 


:/l2  u  2  .  2x1/2 

p  (kx,ky;z)  =  Ae  ~  x  ~  y  ' 


2.14 


At  the  surface  of  plate  A,  the  velocities  of  the  air  and  the  plate  are  equal,  so  that,  from 
Euler’s  equation  with  harmonic  time  dependence,  we  have 


Taking  the  Fourier  transform  yields 


Using  equation  2.10  in  equation  2.18, 


-i4Fp0Qr  irk2k  2  k  2^1/2 

p  (kx,k  ;z)  = - 7 — 2  2  1/2  eH  x  y  ' 

x  y  Mp(k -kx -kyz ) 1/z 

oo  oo  r(2m+l)7tx0-|  r(2n+l)7ty0i 


II 

m=0  n=0 


>1  r 

JcosL 


1^  JC°SL 

l  “ran2  -  0,2  ] 


Imn(kxJCy) 


This  is  the  pressure  radiated  by  plate  A  into  a  free  space,  as  a  function  of  wavenumber, 
kx  and  ky,  at  a  distance  z  from  the  plate.  The  pressure  driving  plate  B  will  be  the 

"blocked"  pressure  which  is  twice  the  free-ficld  pressure  given  in  equation  2.19. 
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2.3.3.  Vibration  Response  of  Plate  B 

The  response  of  plate  B  to  the  pressure  on  its  surface  is  given  by 

oo  o© 

wb(x',y')  =  2  j  j  p(x,y;d)  G(x’,y7x,y)  dx  dy  ,  2.20 

-OO  -oo 

where  wb(x',y')  is  the  velocity  of  plate  B,  p(x,y;d)  is  the  free  field  pressure  at  the 

surface  of  plate  B,  d  is  the  separation  distance  between  the  plates,  and  G(x',y'/x,y)  is 
the  Green's  function  which  satisfies 

DbV4G  -  Pbhbo)2G  «  8(x-x')  8(y-y')  .  2.21 

Equation  2.21  is  the  classical  plate  equation  with  harmonic  time  dependence  and 
pressure  of  unit  amplitude  acting  on  the  plate  at  the  point  x’,y'.  Therefore,  for  the 
infinite  plate,  plate  B,  the  Green's  function,  G(x',y'/x,y)  is  the  velocity  of  the  plate  at 
x',  y'  due  to  a  unit  force  applied  at  a  point  x,y.  Since  the  plate  is  infinite,  the  Green's 
function  depends  only  on  the  separation  between  the  point  source  and  the  receiver,  and 
not  on  each  location.  Equation  2.20  can  be  written  as: 

OO  OO 

wb(x',y')  =  2  j  j  p(x,y;z)  G(x-x'.y-y')  dx  dy  .  2.22 

-OO  -OO 

Equation  2.22  is  the  convolution  integral  and  its  Fourier  transform  is 


Wb(kxJty)  2  p(kxJk.y,z)  G(kx,ky) 


2.23 


For  an  infinite  plate,  it  can  be  assumed  in  equation  2.21  that  x'  =  y'  =  0,  so  that 
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j*t  ^4  ^4  v 

^b(  7  +  2-t-^  +  T-j)0  ■  pbhb0)2O  =  5(x)  5<y)  . 

3x4  3x23y2  ®y 


2.24 


Taking  Fourier  transforms  on  x  and  y,  and  solving  for  the  Green's  function  yields 


^^■Dj^  +  kyV-kp4] 

where  kp^  .the  free  structural  wavenumber,  is  given  by 


2.25 


,  4  Pbhb  2 

k„  = - 0) 

P  Db 


2.26 


Using  equations  2.19  and  2.25  in  equation  2.23,  yields 


-i8FPoe.2  <.i(k2-kx2-ky2)1/2x 

WblkxJty) "  “p  Db[(kx2+ky2)2-kp4](k2-kx2-ky2)1/2 

r(2m+l)icxol  r(2n+l)7ty0-| 

cosL — tz  J^L  c  J 

“  2  ,  ^  ^mn^X’^y^ 

m=0  ri)  t“mn  1  J 


oo  oo 


II 


2.27 


Taking  the  inverse  Fourier  transform,  yields  the  final  solution  for  the  response  of 
plate  B, 


-i8Fp0o>2  Y  V  C°! 


r(2m+l)7tx0-|  r(2n+l)7ty0-| 

L  — e;  JcosL  l;  J 

2  ,.2 


wb<x’y>  * _ La  La 


(2x)2MpDb  “  “  [com„‘-o)‘] 

°°  -  i(k2-kx2-k  2)1/2z .  „  ,  . 

f  f  _J _ JL- _ _  ikx  ikyy  dk  ~ 

L  1  [(k^V^2-^2^  ' 


X  my 


2.28 
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Equation  2.28  has  three  singularities;  one  at  the  resonance  frequencies  of  plate 
A ,  comn  =  co,  one  at  the  acoustic  wavenumber  in  the  air  between  the  plates, 

2  2  2  4  2  2  2 

k  =  kx  -ky  ,  and  one  at  the  free  bending  wavenumber  in  plate  B,  kp  =  (kx  +ky  )  . 

In  order  to  solve  this  equation  these  singularities  must  be  addressed. 

At  resonance  frequency,  the  response  of  plate  A  becomes  infinite  because  there 
is  no  damping  in  the  model.  In  reality,  there  is  always  some  damping,  which  is  added  to 
the  model  by  making  the  bending  rigidity  of  plate  A  complex, 

Da*=  Da(l+illa)  ,  2.29 

where  qa  is  the  damping  loss  factor  of  plate  A. 

The  singularity  at  the  acoustic  wavenumber  results  from  taking  a  Fourier 
transform  over  the  infinite  spatial  domain  of  plate  B.  With  no  dissipation,  the  acoustic 
pressure  waves  in  the  air  will  propagate  parallel  to  plate  B  out  to  infinity,  so  that  energy 
will  sum  to  infinity  when  a  sum  is  taken  over  all  space  in  the  infinite  Fourier  transform. 
In  reality,  there  will  be  some  dissipation  of  the  acoustic  energy  propagating  parallel  to 
plate  B.  Also,  energy  that  propagates  parallel  to  plate  B  should  couple  less  to  the 
response  of  plate  B  then  energy  that  propagates  directly  across  the  space  between  the 
plates  and  is  incident  on  plate  B  at  non-grazing  angles.  Therefore  dissipation  can  be 
added  to  the  acoustic  propagation  to  make  the  model  more  realistic  and  remove  the 
singularity  at  the  acoustic  wavenumber,  without  impacting  the  results  at  non-acoustic 
wavenumbers. 

Pierce  (5  eqn.  10-8.9b)  gives  the  acoustic  wavenumber  with  absorption  as: 

*  co 

k  =  c  +  icXgj  , 


2.30 


where  the  classical  absorption  constant  is  given  by  (from  Pierce  (5  eqn.10-2.12)) 


a)2 

acl'  = 


2p0 


V  3  +  ,,  +  Pr  / 


2.31 


where  p.  is  the  viscosity  of  air, 
pc  is  the  density  of  air, 
pB  is  the  bulk  viscosity  of  air, 
y  is  the  ratio  of  specific  heats  =  1.4  for  air,  and 
Pr  is  the  Prandtl  number  for  air. 

Ignoring  effects  of  bulk  viscosity  (pB)  and  using  equation  2.31  in  equation  2.30, 


Equation  2.32  was  used  for  the  acoustic  wavenumber  in  equation  2.28. 

The  singularity  at  the  free  bending  wavenumber  in  plate  B  can  be  removed  by 
adding  damping  to  the  plate  in  a  manner  similar  to  the  damping  added  to  plate  A.  This 
damping  has  the  same  effects  as  the  dissipation  added  to  the  air  between  the  plates.  With 
damping,  the  bending  wavenumber  can  be  written  as 

(kp*)4-kp4(l+ir|b)  .  2.33 


where  Tjj,  is  the  damping  loss  factor  of  plate  B. 


2.3.4.  The  Discrete  Inverse  Fourier  Transform 
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To  evaluate  the  integral  in  equation  2.28  using  numerical  methods,  a  two 
dimensional  discrete  inverse  Fourier  transform  was  developed  based  on  the  continuous 
infinite  inverse  Fourier  transform  as  defined  by  Brigham  (6), 


oo  oo 


f(x,y)  =  J  J  F(kx,ky)  e  lkxx  e  ikyy  dkx  dk  . 


-oo  -oo 


This  integral  is  approximated  by  limiting  the  range  of  integration,  so  that 


Kx/2  KJ2 

f(x,y)  -  J  /  F(kx,ky)  e  ^x  e  ik 


-Kx/2  -Ky/2 


yy  dkx  dky 


where  Kx  and  Ky  are  at  least  twice  the  value  of  kp  as  calculated  from  equation  2.26  or 

twice  the  acoustic  wavenumber  k  which  ever  value  is  larger  at  the  frequency  of  interest. 

These  limits  were  selected  so  that  the  range  of  integration  includes  all  of  the  points 

where  singularities  would  exist  in  the  absence  of  added  damping  and  dissipation  and  so 

i(k  -k  ^-k 

that  the  exponential  term  e  x  y  decays  to  a  small  number  at  the  limits  of 
integration.  Equation  2.3S  can  be  written 


Kx  Ry 

f(x,y)  =  e  -lxKx/2  e  'iyKy/2  J  J  F(kx  -  Kx/2,ky  -  Ky/2)  e  ^x*  e  ikyy  dkx  dky  . 

0  0 


Kx2it 


This  integral  can  be  evaluated  using  the  rectangular  rule  with  spacing  h^  =  and 
Ky27t 

hky  =  N  so 
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f(X'y)  *  (X?  '  '  iy  Ky/2hkxhky 

N-l  N-l 

X  X  c  lxm^kx  c  iynhky  F(  mh^-  Kx/2,nhky  -  Ky/2)  .  2.37 

m=0  n=0 

2jq'x  2ftjv 

By  setting  x  =  ~k~  and  y  =  “j^  where  jx  =  0,1,...N-1  and  jy  =  0,1.— N-l.  Equation 
x  y  j 

2.37  becomes 


27tjx2rcjv 
f  KX  *  Ky 


)= 


(-l)jx(.i)jyhkxhky 


N-l  N-l  K27t)2jxm  i(2fl)2jyn 

X  Xc  N  «  N  F(mhkx-Kx/2,nhky-KJ/2)  . 


1  ♦ 

Thus,  after  scaling  by  the  factors  "^2  (”1)  (-1)  ^  hkxhky>  the  discrete  inverse 


Fourier  transform  given  by  equation  2.38  approximates  the  continuous  inverse  Fourier 
transform  given  by  equation  2.35. 

To  calculate  the  plate  response,  the  infinite  summation  in  equation  2.28  can  be 
solved  by  starting  the  summation  at  the  resonance  frequencies,  (0^,  which  are  nearest 

the  frequency  of  interest  With  this  approach,  the  summation  converges  more  rapidly 
than  beginning  with  the  lowest  mode  numbers.  This  is  seen  by  observing  equation  2.28. 

When  a  resonance  frequency  is  chosen  which  is  not  near  the  resonance  frequency  of 
interest  the  (co^2  -  to2)  term  in  the  denominator  becomes  large  and  the  contribution  to 


the  sum  becomes  small.  When  solving  equation  2.28  an  iterative  process  was  used.  The 


22 


summation  was  taken  using  the  resonance  frequencies  near  the  frequency  of  interest 
until  the  contribution  to  the  summation  became  negligible.  The  same  process  can  be  used 
to  solve  the  infinite  summation  in  equation  2.9. 

2.3.5.  Transmission  Loss  between  Plate  A  and  Plate  B 


Transmission  loss  is  given  by 


TL  —  10  log  jq 


(V 

(Ab? 


2.39 


where  Aa  is  the  average  of  the  magnitude  of  the  velocity  of  plate  A  vibration,  computed 
from  equation  2.9  and  Afi  is  the  average  of  the  magnitude  of  the  velocity  of  plate  B 
vibration,  computed  from  equation  2.28. 

For  purposes  of  computation,  the  transmission  loss  for  the  vibration  response  of 
plate  B  is  computed  in  the  shadow  zone  of  plate  A  as  illustrated  in  Figure  2.1.  By 
computing  the  vibration  amplitude  at  N  points  on  plate  A  and  M  points  on  plate  B,  the 
transmission  loss  is  determined  using  the  following  equation: 


1  N 

b  X<aa‘>2 

TL  =  10  log10  — “ 

Si  X(abj)2 

j=l 


N  and  M  are  chosen  so  that  there  are  three  points  per  wavelength  at  the  highest 
frequency. 
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Chapter  3 
RESULTS 

The  equations  developed  in  the  previous  chapter  are  used  to  predict  airborne 
transmission  losses  for  plate  configurations  which  are  representative  of  typical  machine/ 
foundation  arrangements.  Table  3.1  shows  the  input  values  used  in  the  equations.  For 
various  input  parameters,  surface  plots  of  the  displacements  of  plates  A  and  plate  B  were 
obtained  to  show  the  motion  of  the  plates  and  graphs  of  airborne  transmission  loss 
between  plates  A  and  B  as  a  function  of  frequency  are  presented. 

3.1.  Surface  Plots 

Figures  3.1  through  3.6  are  surface  plots  of  the  motion  of  plates  A  and  B  for 
various  frequencies  and  damping.  These  surface  plots  are  essentially  a  snap  shot  of  the 
displacement  of  the  plate.  The  surface  plots  for  plate  A  were  obtained  using  equation  2.9 
in  a  FORTRAN  computer  program.  The  surface  plots  for  plate  B  were  obtained  using 
equation  2.28  in  a  FORTRAN  computer  program.  The  discrete  inverse  Fourier 
transform  as  described  in  equation  2.38  was  used  to  solve  the  integral  part  of  equation 
2.28. 

In  Figures  3.1  through  3.6  the  figure  with  the  "a"  suffix  is  the  surface  plot  of 
plate  A  and  the  figure  with  the  "b"suffix  is  the  surface  plot  of  plate  B.  The  x  and  y 
ordinates  are  the  dimensions  of  the  plate  in  meters  (m).  The  vertical  axis  is  the  amplitude 
of  displacement  of  the  plate  in  meters.  The  parameters  used  to  make  the  predictions 
presented  in  Figures  3.1  through  3.6  are  listed  in  Table  3.1  and  shown  on  the  surface 
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Table  3.1 

Input  parameters  for  calculation  of  surface  plots  and  transmission  loss  curves 


Parameters  of  Finite  Plate  (Plate  A) 

dimensions  of  plate . 

..1^  =  1.0  m  Ly  =  1.0  m 

driving  point . 

.x0  =  0.0  m  y0  =  0.0  m 

driving  force . 

..F  =  1.0  N 

plate  thickness . 

,.ha  =  0.001587  m 

plate  density . 

•Pa  =  7860.0  kg/m2 

modulus  of  elasticity . . 

...Ea  =  200  x  109  N/m2 

bending  rigidity . 

.Da  =  72.73  N/m 

damping  coefficient . 

...ria  =  0.001 

Parameters  of  Infinite  Plate  (Plate  B) 

plate  thickness . 

..hb  =  0.00158 7  m 

plate  density . 

..pb  =  7860.0  kg/m2 

modulus  of  elasticity . 

...Eb  =  200  x  109  N/m2 

bending  rigidity . 

,.Db  =  72.73  N/m 

damping  coefficient . 

...T)b  =  0.001 

Parameters  for  Fluid  (Air) 

density . 

..pb  =  1.204  kg/m2 

viscosity . 

...p  =  0.0000184  kg/ms 

speed  of  sound . 

,...c  =  343.0  m/s 

Prandtl  number . 

...Pr  =  0.706 

specific  heat  ratio . 

....y  =  1.4 

distance  between  plates... 

. z  =  0.1  m 
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plots.  Surface  plots  were  generated  for  a  low  resonant  frequency,  a  low  non-resonant 
frequency  and  a  high  resonant  frequency,  for  both  high  and  low  damping.  Plate  A  is  a 
square  plate  with  a  sinusoidal  force  of  1  Newton  (N)  (0.2248  lb.)  applied  in  the  center. 
Both  plates  are  steel.  The  edges  of  plate  A  are  at  zero  amplitude  as  required  by  the 
simple  supports  at  the  edges. 

Figures  3.1,  3.2,  3.3,  and  3.4  show  plate  A  and  plate  B  excited  at  37.93  Hz 

which  from  equation  2.7  is  the  n  =  0,  m  =  1  resonant  mode  of  plate  A.  In  figures  3.1 
and  3.2  both  plates  are  highly  damped  (qa  =  =  0.1)  and  in  Figure  3.3  and  3.4  both 

plates  are  lightly  damped  (na  =  =  0.001).  Comparing  Figures  3.1  and  3.3  it  can  be 

seen  that  the  amplitude  of  the  undamped  plate  A  is  much  larger  than  the  highly  damped 
plate  A.  This  illustrates  the  well  known  effect  of  damping  at  reducing  the  resonant 
response  of  structures.  Comparing  the  responses  of  the  damped  and  undamped  plate  A 
(Figures  3.1  and  3.3),  it  can  be  seen  that  the  mode  shapes  arc  the  same  even  though  the 
maximum  amplitude  of  the  plate  vibration  is  different 

The  symmetry  of  the  waves  excited  in  plate  B  and  propagating  outward  from  the 
shadow  zone  of  plate  A  can  be  seen  clearly  in  Figure  3.2.  In  this  figure,  the  exponential 
decay  due  to  damping  can  also  be  seen.  The  amplitude  of  the  waves  in  the  middle  of 
plate  B  directly  beneath  the  top  plate  is  large  and  the  waves  decay  exponentially  as  they 
spread  towards  the  edge  of  the  plate.  In  Figure  3.4  for  the  lightly  damped  plate  B  it  can 
be  seen  that  the  plate  is  fully  excited.  Also,  because  of  the  decrease  in  the  amplitude  of 
vibration  of  plate  A,  the  vibration  of  plate  B  in  the  shadow  zone  is  lower  with  high 
damping  than  with  low  damping. 

Figures  3.5,  3.6,  3.7,  and  3.8  show  plate  A  and  plate  B  excited  at  1 10  Hz, 
which  is  not  a  resonance  frequency  of  plate  A.  It  is  between  the  m  =  2,  n  =  0  resonance 
at  98  Hz  and  the  m  =  2,  n  =  1  resonance  at  129  Hz.  The  mode  shape  of  the  damped  and 
undamped  plate  A  (Figures  3.5  and  3.7)  are  similar  to  each  other.  The 
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AMPLITUDE 


Figure  3.2  Plate  B,  low  resonance, high  damping  (i 
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A,  low  resonance,  low  damping  (f=37.9  Hz  and  t\  =  0.001) 


■258.5 


Figure  3.4  Plate  B,  low  resonance,  low  damping  (f=37.9  Hz  and  B  =  0.001) 
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Figure  3.6  Plate  B,  low  non-resonance,  high  damping  (f=l 
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Figure  3.7  Plate  A,  low 
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difference  between  the  amplitudes  of  the  damped  and  undamped  plate  A  is  not  as  great 
as  at  the  resonance  frequency  (Figures  3.1  and  3.3).  This  is  because  damping  in  the 
plate  is  less  effective  at  reducing  structural  response  of  a  structure  off  resonance  than  at 
resonance. 

The  radial  symmetry  of  the  waves  excited  in  plate  B  can  again  be  seen  in  Figure 
3.6  with  the  high  damping.  The  exponential  decay  of  the  waves  toward  the  edges  of 
plate  B  produced  by  the  damping  can  also  be  seen.  The  response  in  the  shadow  zone 
shown  in  Figure  3.6  is  negative,  opposite  that  shown  in  Figure  3.1.  Again  with  light 
damping,  plate  B  is  excited  both  in  and  out  of  the  shadow  zone,  as  shown  in  Figure  3.8. 

Figure  3.9  is  a  surface  plot  of  the  inverse  Fourier  transform  of  equation  2.19 
using  the  same  input  values  as  for  Figure  3.1  and  3.2.  This  is  the  pressure  incident  on 
plate  B  at  a  low  resonance  frequency  and  high  damping  (37.9  Hz  and  q=0.01).  It  can  be 
seen  in  this  figure  that  the  pressure  distribution  on  plate  B  closely  follows  the  mode 
shape  of  plate  A  as  shown  in  Figure  3.1.  Figure  3.10  is  a  surface  plot  of  the  inverse 
Fourier  transform  of  equation  2.25  using  the  same  inputs.  This  is  the  Green's  function 
of  Plate  B.  Equation  2.23  is  the  response  of  plate  B  and  is  obtained  by  combining 
equations  2.19  and  2.25  and  multiplying  by  2.  Similarly,  if  each  point  on  Figures  3.9 
and  3.10  are  combined  and  multiplyed  by  2,  the  result  is  Figure  3.2,  the  response  of 
plate  B.  Figures  3.11  and  3.12  are  the  pressure  and  Green's  functions  respectively  for 
1 10.0  Hz  and  q=0.1.  These  two  figures  can  be  combined  in  the  same  manner  to  obtain 
the  plate  response  shown  on  Figure  3.6.  Comparing  the  pressures  shown  in  Figures  3.9 
and  3.1 1,  the  effects  of  the  resonant  response  of  the  plate  on  the  increase  in  the  pressure 
incident  on  plate  B  in  the  shadow  zone  can  be  seen.  At  resonance,  the  incident  pressure 
is  higher.  Figures  3.10  and  3.12  show  the  effects  of  free  propagation  in  the  infinite  plate 
outside  the  shadow  zone.  At  the  higher  frequency,  the  wavelength  for  free  propagation 
is  smaller. 
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Figure  3.9  Pressure  incident  on  plate  B,  low  resonance,  high  damping  (f=37.9Hz  and  r|  =  0. 1 ) 
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Figure  3.10  Green's  function  for  plate  B,  low  resonance,  high  damping  (f=37.9Hz  and  "H  =0.1) 
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Figure  3.1 1  Pressure  incident  on  plate  B,  low  non-resonance,  high  damping  (f=l  10.01  Iz  and  B  =  0.1) 
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Figure  3.12  Green's  function  for  plate  B,  low  non-Tesonance,  high  damping  (f=l  10.0Hz  and  B  =  0.1) 


39 


3.2.  Transmission  Loss  between  plates  A  and  B 

Predictions  of  the  motions  of  plate  A  and  plate  B  obtained  by  using  equations 
2.9  and  2.28  in  equation  2.40  to  predict  the  transmission  loss  in  decibels  from  10  Hz  to 
5000  Hz  are  presented  in  Figures  3.13  through  3.18.  Unless  otherwise  noted,  the 
parameters  listed  in  Table  3.1  are  used,  with  the  separation  distance  between  the  plates, 
the  thickness  of  plate  B,  the  size  of  plate  A,  and  the  damping  varied.  The  values  in  Table 
3.1  were  chosen  to  simulate  a  typical  sub-base  and  foundation  for  a  medium  sized  piece 
of  equipment.  Also  shown  in  Figures  3.13  through  3.18  is  the  theoretical  transmission 
loss  for  simple  resilient  mounts  with  a  resonance  frequency  of  5  Hz  as  computed  in 
equation  2.1  and  shown  in  Figure  2.2.  The  parameters  given  by  Snowdon  (2)  for 
rubber  mounts  were  used  in  equation  2.1. 

In  Figure  3.13  the  effect  on  airborne  transmission  loss  of  changing  the 
separation  distance  between  the  plates  from  0.01  m  to  0.5  meters  is  presented.  Above 
100  Hz  all  of  the  curves,  except  for  the  curve  for  0.5  m,  fall  below  the  transmission  loss 
curve  for  the  mounts.  This  implies  that  the  airborne  path  would  degrade  the  performance 
of  the  isolation  system  even  for  the  largest  separation  distances.  Increasing  the 
separation  distance  increases  the  airborne  transmission  losses  at  frequencies  between 
100  and  IK  Hz.  Below  100  Hz,  the  separation  distance  has  little  impact  on  the  airborne 
transmission  losses,  other  than  at  the  frequencies  where  maxima  and  minima  occur  in 
the  transmission  loss.  The  minima  in  the  transmission  loss  curve  below  100  Hz  may  be 
due  to  a  resonance  associated  with  the  effective  masses  of  the  plates  and  the  effective 
stiffness  of  the  air  between  the  plates.  For  separation  distances  greater  than  0.05  m, 
plate  B  moves  out  of  the  acoustic  nearfield  of  plate  A  at  frequencies  above  IK  Hz,  so 
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Figure  3. 13  Transmission  loss  curves,  changing  separation  distance  between  the  plates 
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that  the  transmission  losses  above  IK  Hz  become  less  sensitive  to  separation  distance 
than  at  frequencies  below  IK  Hz. 

In  Figure  3.14  the  effect  on  airborne  transmission  of  changing  the  thickness  of 
plate  B  from  0.1587  cm  (1/16  in.)  to  1.27  cm  (1/2  in.)  is  shown.  Below  IK  Hz, 
increasing  the  thickness  of  the  receiver  plate  increases  the  airborne  transmission  losses  at 
all  frequencies.  Above  100  Hz,  the  transmission  losses  for  the  0.0157  cm  (1/16  in.) 
thick  plate  are  below  the  transmission  losses  for  the  mounts.  Below  IK  Hz,  the 
transmission  losses  for  the  two  thicker  plates  are  above  the  losses  for  the  mounts. 
Above  IK  Hz,  the  transmission  losses  for  all  three  plates  are  below  the  losses  for  the 
mounts.  Increasing  the  stiffness  of  the  receiving  plate  by  increasing  the  thickness 
increases  the  airborne  transmission  loss  so  that  the  airborne  path  is  not  significant  at 
lower  frequencies  (below  IK  Hz);  however,  at  higher  frequencies  (above  IK  Hz)  the 
airborne  path  remains  a  problem  even  when  the  stiffness  of  the  structure  is  high. 

In  Figure  3.15  the  size  of  plate  A  is  changed  from  1  m  x  1  m  to  0.1  m  x  0.1  m. 
Above  100  Hz,  the  transmission  losses  for  the  two  larger  plates  are  below  the 
transmission  losses  for  the  mounts  and  above  300  Hz,  the  airborne  transmission  losses 
for  the  smallest  plate  is  also  below  the  losses  for  the  mounts.  Decreasing  the  area  of 
plate  A  without  changing  the  separation  distance  increases  the  ratio  of  the  area  on  the 
sides  between  the  plates  available  for  air  to  escape  to  the  surface  area  of  plate  A  and  the 
surface  of  plate  B  in  the  shadow  zone.  Thus,  with  the  smaller  plate,  a  larger  percentage 
of  the  energy  in  between  the  plates  is  permitted  to  escape  which  increases  the  airborne 
transmission  losses  between  the  plates  at  low  frequencies.  At  the  higher  frequencies,  air 
entrapment  becomes  less  of  a  factor,  decreasing  the  effect  of  the  size  of  plate  A  on  the 
airborne  transmission  loss.  Also  the  dip  and  peak  in  the  airborne  transmission  loss  curve 
for  the  large  (1  m  x  1  m)  plate  are  reduced  when  the  size  of  plate  A  is  decreased.  This 
may  be  a  damping  effect  produced  by  an  increase  in  the  percentage  of  the  energy 
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Figure  3.14  Transmission  loss  curves,  changing  the  thickness  of  plate 
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Figure  3.15  Transmission  loss  curves,  changing  the  size  of  plate  A 
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between  the  plates  lost  to  propagation  out  from  between  the  plates  through  the  openings 


at  the  sides  of  the  plates. 

In  Figures  3.16  and  3.17  the  damping  in  plates  A  and  B  was  changed  for  two 
different  separation  distances  between  the  plates.  In  Figure  3.16  the  distance  between 
the  plates  is  0.1  m  and  the  damping  is  increased  from  0.001  (a  typical  number  for  an 
undamped  steel  plate)  to  0.1  (a  typical  number  for  steel  with  damping  applied).  In 
Figure  3.17  the  distance  between  the  plates  is  0.01  m  and  the  damping  is  varied  in  the 
same  manner.  For  both  of  these  graphs  the  airborne  transmission  loss  curve  for  the 
highly  damped  plate  follows  the  same  trend  as  for  the  lighdy  damped  plate;  however, 
with  damping,  the  curves  are  much  flatter.  This  shows  the  effect  of  damping  which 
works  primarily  at  resonances.  For  both  separation  distances,  damping  increases  the 
airborne  transmission  losses  at  63  Hz  where  the  airborne  transmission  loss  with  light 
damping  is  below  the  transmission  loss  curve  for  the  mounts.  However,  above  100  Hz, 
where  the  airborne  transmission  loss  curve  does  not  display  the  effects  of  resonances, 
damping  has  less  effect  on  the  airborne  transmission  losses. 

Adding  damping  to  only  plate  B  produced  the  results  shown  in  Figure  3.18  for  a 
separation  distance  of  0.1  m.  The  damping  in  plate  B  has  an  effect  only  at  the 
frequencies  below  100  Hz  where  a  peak  and  a  dip  in  the  airborne  transmission  loss 
curve  occurs.  Again  the  damping  appears  to  damp  the  resonant  effect  that  occurs 
between  the  plates  and  the  air  cushion  between  the  plates. 
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Figure  3.16  Transmission  loss  curves,  changing  the  damping  in  plates  A  and  B 
for  a  plate  separation  distance  of  0. 1  m 
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Figure  3.17  Transmission  loss  curves,  changing  the  damping  in  plates  A  and  B 
for  a  plate  separation  distance  of  0.01  m 
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Figure  3.18  Transmission  loss  curves,  changing  the  damping  in  plate  B 
for  a  plate  separation  distance  of  0. 1  m 


Chapter  4 
CONCLUSIONS 
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An  analytic  model  for  the  airborne  transmission  loss  between  a  vibrating  finite 
plate  and  a  parallel  infinite  plate  was  developed.  Comparison  of  the  predicted  airborne 
transmission  losses  with  predicted  transmission  losses  for  resilient  mounts  revealed  that 
airborne  transmission  losses  were  often  less  than  the  losses  through  the  mounts  at 
frequencies  above  100  Hz,  implying  that  the  airborne  transmission  path  may  be 
significant  in  many  well-designed  single-stage  resilient  mounting  systems.  Below 
100  Hz,  the  transmission  losses  through  the  mounts  are  below  the  airborne  transmission 
losses,  except  at  frequencies  where  resonances  between  the  plates  occur. 

Increasing  the  separation  distance  between  the  plates,  increasing  the  thickness  of 
the  receiver  plate  and  decreasing  the  size  of  the  source  plate  increased  the  airborne 
transmission  losses  between  the  plates.  Increasing  the  damping  in  the  plates  was  most 
effective  in  reducing  peaks  and  dips  in  the  airborne  transmission  loss  curves  below 
100  Hz. 

The  basic  assumptions  used  in  developing  the  analytic  model  are  that  the 
impedances  of  the  plates  are  much  larger  than  the  impedance  of  the  air  between  the 
plates,  there  is  no  feedback  from  the  receiver  plate  to  the  source  plate  and  that  the 
pressure  generated  in  the  plane,  but  outside  the  area,  of  the  source  plate  has  a  negligible 
effect  on  the  airborne  transmission  losses.  In  future  research,  these  assumptions  should 
be  relaxed,  particularly  the  feedback  assumption  where  standing  waves  between  the 
plates  should  decrease  the  airborne  transmission  losses.  Also,  relaxing  the  small  air 
impedance  assumption  would  increase  the  coupling  between  the  air  and  the  plates 
increasing  the  vibratory  response  of  the  plates  and  decreasing  the  airborne  transmission 
losses.  Experimental  verification  of  the  analytic  model  is  needed  to  certify  the  conclusion 
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presented  above.  Also,  experimental  and  theoretical  analysis  of  treatments  to  increase  the 
airborne  noise  transmission  losses  between  structures,  such  as  decoupling  coating  and 
damping,  is  needed. 
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Appendix  A 
CALCULATION  OF 


In  this  appendix,  the  integral 
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Wkx’ky)  =  J  /  cos [(2X)KX] cos [(2n^)7Cy]  c  c  'V  dx  dy  A.l 


as  given  in  equation  2.1 1  is  evaluated. 

Since  cos  [  ]  and  cos  ["  ]  are  even  and  defined  only  for 
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which,  when  integrated,  yields 
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COMPUTER  MODEL 
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In  order  to  solve  the  equation  of  motion  of  both  plate  A  and  plate  B  and  to 
compute  the  transmission  loss  between  the  two  plates,  three  separate  computer  programs 
were  developed. 

Program  PRA  (Plate  Response  of  plate  A)  computes  the  mode  shape  of  plate  A. 
Program  PR  (Plate  Response)  computes  the  response  of  plate  B.  Program  PRTL  (Plate 
Response  -  Transmission  Loss)  computes  the  transmission  loss  between  plate  A  and 
plate  B.  PRTL  is  broken  down  into  three  parts  for  three  different  frequency  ranges: 
PRTL100  for  the  frequency  range  from  10  to  100  Hz,  PRTL1000  for  the  frequency 
range  from  100  to  1000  Hz  and  PRTL5K  for  the  frequency  range  from  1000  to  5000 
Hz.  The  source  code  for  each  of  these  programs  is  in  Appendix  C. 

The  inputs  to  all  of  these  programs  are  taken  from  the  same  input  file  format.  A 
sample  input  and  its  corresponding  input  variables  are  included  in  Appendix  C.  Not  all 
of  these  inputs  are  used  in  each  of  the  programs.  See  the  corresponding  write  up  below 
for  which  variables  are  needed  for  each  program. 

B.l.  Program  PRA 

This  program  solves  equation  2.9  to  obtain  the  mode  shape  of  plate  A.  The 
inputs  to  this  program  are: 

Lx  and  Ly  Dimensions  of  plate  A 

Xo  and  Yo  The  point  where  the  input  force  is  applied  to  plate  A 

ROA  Density  of  plate  A 

ACHA  Thickness  of  plate  A 

ETAA  Damping  coefficient  for  plate  A 

EA  Modulus  of  Elasticity  of  plate  A 
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The  output  from  this  program  is  a  surface  plot  of  the  mode  shape  of  plate  A.  The 
plotting  program  used  to  obtain  these  plots  is  a  subroutine  called  ARL_HIDE  and  is 
available  on  the  ARL  Penn  State  VAX  computer. 


B.2.  Program  PR 


This  program  uses  equation  2.28  to  obtain  the  vibration  response  of  plate  B  to 
excitation  by  airborne  transmission  from  plate  A.  The  integral  portion  of  equation  2.28 
is  solved  using  a  2-dimensional  discrete  inverse  fast  Fourier  transform  program  called 
F2T2B  in  the  IMSL  subroutine  library,  also  available  on  the  ARL  Penn  State  VAX 
computer.  In  equation  2.28  this  is  an  infinite  2-dimensional  integral;  however,  in  order 
to  solve  this  integral  using  the  computer  the  finite  integral  is  approximated  by  a  finite 
series  using  equation  2.38. 

The  choice  of  wavenumber  range  (KX1  and  KYI)  and  the  sample  rate 
determines  the  size  of  plate  B  by  the  following  equations. 


2N  2N 

x  =  KXl  y  =  KYI 


B.l 


where  x  and  y  are  the  dimensions  of  plate  B. 

The  wavenumber  range  must  be  chosen  judiciously  with  respect  to  the  frequency 
of  interest  because  of  the  relationship  of  the  wavenumber  and  frequency  in  equation 
2.28.  At  lower  frequencies  a  lower  wavenumber  range  can  be  used,  therefore  a  larger 
area  of  plate  B  can  be  observed  using  equation  B.l.  At  higher  frequencies  a  higher 
wavenumber  range  must  be  used  and  a  smaller  area  of  plate  B  can  be  observed  using  the 
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same  sampling  rate  (N).  A  higher  sampling  rate  will  solve  this  problem  however,  the 
computer  time  also  increases. 

The  inputs  to  this  program  are: 

Lx  and  Ly  Dimensions  of  plate  A 

Xo  and  Yo  The  point  where  the  input  force  is  applied  to  plate  A 

ROA  Density  of  plate  A 

ROB  Density  of  plate  B 

ACHA  Thickness  of  plate  A 

ACHB  Thickness  of  plate  B 

ETAA  Damping  coefficient  for  plate  A 

ETAB  Damping  coefficient  for  plate  B 

EA  Modulus  of  Elasticity  of  plate  A 

EB  Modulus  of  Elasticity  of  plate  B 

Z  Distance  between  the  plates 

PRA  Prandel  number  for  the  fluid 

GAMMA  Specific  heat  ratio  for  the  fluid 

MU  Viscosity  of  the  fluid 

NU  Poison's  ratio 

C  Speed  of  sound  in  the  fluid 

RO  Density  of  the  fluid 

KX1  Upper  wavenumber  limit  in  the  x  direction 

KYI  Upper  wavenumber  limit  in  the  y  direction 

The  output  of  this  program  is  a  surface  plot  of  the  vibration  response  of  plate  B. 
ARL.HIDE  is  also  used  to  obtain  these  plots. 


B.3.  Program  PRTL 
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Program  PRTL  uses  the  same  algorithms  used  in  programs  PRA  and  PR  to 
obtain  the  motions  of  plates  A  and  B.  Equation  2.40  is  then  used  to  obtain  the 
transmission  loss  between  the  two  plates  in  the  shadow  zone  of  plate  A. 

This  program  has  been  broken  down  into  three  parts;  PRTL  100,  PRTL 1000  and 
PRTL5K  in  order  to  change  the  sampling  rate  and  wavenumber  range  for  the  frequency 
range  covered  by  each  PRTL  program. 

The  input  to  these  programs  are  the  same  as  for  program  PR. 

The  output  from  these  programs  is  a  table  of  transmission  losses  as  a  function  of 
frequency. 


Appendix  C 

FORTRAN  COMPUTER  PROGRAMS 


PROGRAM  PKA 

C  PLATE  RESPONSE  OF  PLATE-A,  FINITE  PLATE  _  _  _ _ 

C  WRITTEN  BY  MICHAEL  F  SHAW 

INTEGER’  OMEGA  I OMEGAS,  OMEGAtT  INC,  "emTeN,  N,  I,  J, 

AN1.N2.N11  .  .  .  _ 

C  NSHOULD  BE  A  PRODUCT  OF  SMALL  PRIMES 

PARAMETER  (N-12B)  .  .  _ 

COMPLEX  DSTAR , OMEGAMNSO . B ( N , N )  ,BT(N,N)  ,WA(N,N)  ,WA2(N,N) , 

_  ABT2JN.N)  ...  _ _ 

REAL  VMAX,VMIN,VINC,WAREAL(N+1,N+1)  ,WORK1{4000) ,W0RK2(100) 

_ REAL  DA ,  ETAA .  ETAB , LX , LY ,  ROA ,  ROB .  AG!  A .  AQ1B ,  WA2REAL  ( N+l ,  N+l )  _ _ 

AKXl, KXl.F.MP.XO.yO.NU, PI, OMEGA, GAMMA.RO, MU, C.EA.EB, 2, 

_ A0MEGA1  .QMEGA2  .QMEGAMNO ,  QMEGAMNl ,  QMEG  AMN2  . . 

INTRINSIC  CMPLX 


C  VARIABLE  LIST 

C  _  .  .  . 

C  DA  -  BENDING  RIGIDITY  OF  THE  FINITE  PLATE 

C._.  ETAA  -  DAMPING  COEFICIENT  FOR  THE  FINITE  PLATE  _ 

C  LX  4  LY  -  DIMENSIONS  OF  THE  FINITE  PLATE 

C_  ROA  -  DENSITY  OF  THE  FINITE  PLATE  .  ... _ 

C  ACHA  -  THICKNESS  OF  WE  FINITE  PLATE 

C . RQ  t  DENSITY  OF  FLUID  _  _ _ 

C  XO  4  YO  -  TOE  POINT  WHERE  TOE  INPUT  FORCE  F  IS  APPLIED 

C__  NU  -  POISON'S  RATIO  . . 

C  EA  -  SHEER  MODULUS  OF  TOE  FINITE  PLATE 


PARAMETER  (F-1.0) 

C_  INPUT  FORCE  TO  TOE  FINITE  PLATE  IS  ASSUMED  TO  BE  A  UNIT  FORCE  (1) 
OPEN  (UNIT-25,  FILE- *  PR . DAT ' ,  STATUS- ‘ OLD * ) 

_ _  OPEN  (UNIT-26,  FILE-'PRA.OUT' ,  STATUS-'NEW' )  _ _ 

C  READ  INPUT  DATA  FROM  PR. DAT 

_  READ  (25,15)  LX.LY.XO, YO.ACHA.ACHB  _  _ 

15  FORMAT( 4F8 . 3 , 2F8 . 6 ) 

.  .  _  READ  (25.25)  ROA, ROB, RO. ETAA. ETAB, GAMMA  ..  _  .  _  ... 

25  FORMAT( 2FB . 2 , F8 . 3 , 2F8 . 7 , F8 . 6 ) 

READ  (25,35)  MU.C.EA.EB.KXl.KYl  .  _  _ 

35  FORMAT( F8 . 7 , FB . 3 , 2E8 . 3 , 2F8 . 2 ) 

READ  (25,45)  NU, 0MEGA1, QMEG A2, INC, 2, OMEGA  .  .  _ 

45  P0RMAT(F8. 4 , F8. 1 , F16 . 1 , I8,F8 . 3 ,F12 . 3) 

C  COMPUTE  THE  VALUE  OF  PI  ....  _  ...  . 

PI-2. 0*ASIN(  1.0) 

C  CQMPUT  TOE  BENDING  RIGIDITY  _  .  .  ........  _  ...... 

DA-(EA*AQIA**3. 0)/( 12 . 0« ( 1 . 0-NU»*2 . 0) ) 

DSTAR -CMP  LX  ( DA ,  DA*  ETAA ) 

MP-ROA*ACHA*LX*LY 

WRITE(26,11)  LX.LY  .._... . 

11  FORMAT ( '  PARAMETERS  OF  FINITE  PLATE', /,'  DIMENSIONS  OF  PLATE' 

_ 4,'  (LX)  ' ,F8. 3. '  (LY)  ' ,F8 . 3 , '  METERS')  ...  .  . 

WRITE(26,21)  XO,YO 

21  FORMATC  DRIVING  POINT  XO  ',FB.3,'  YO  ',FB.3,' 

4  METERS') 

WRITE(26.31)  F 

31  FORMATC  DRIVING  FORCE  ',F8.3,'  N') 

WRITE ( 2  6 , 4 1 ) ACHA 

41  FORMATC  PLATE  THICKNESS  \F8.6,'  M') 

WRIT£(26,51)  ROA  .  _. 

51  FORMATC  PLATE  DENSITY  ',F8.3,'  KG/CU-METERS ' ) 

WRITE (26, 61 )  CA 

61  FORMATC  MODULUS  OF  ELASTICITY  ',E0.3,'  N/SQ-METERS ' ) 

_  WRITE(26.71)  DA  ....._  _ 

71  FORMATC  BENDING  RIGIDITY  ‘.E12.4) 

WRITE! 26, 81)  ETAA  _  _  _ _ 

81  FORMATC  DAMPING  COEFFICIENT' ,F8. 4) 

C  TOE  FREQUENCY  RANGE  IS  TAKEN  FROM  OMEGA 1  TO  OMEGA2  IN  INCREMENTS 


CORRESPONDING  TO  INC 


C  OMEGAB-10*IFIX(  ALOG10(OM£GAl ) ) 

C _ OMEGAT-INC*10*IEIX(ALOG10(OMEGA2)) _ 

C  DO  1000  OMEGAI-OMEGAB.OMEGAT 

C _ QMEGA-2.0*PI? (IQ.  Of? (QMEGAI/(  INC*10. 0)  )  ) _ _ _ 

DO  100  I-l.N 

_ DO  100  Jrl.N _  _ _ _ 

BT(I,J)-0.0 

_ WAREAL(I,J)=Q.Q _ _ _ 

100  CONTINUE 

_ WHITE (26 ,12)  OMEGA/ (2 . 0*PI ) _ _ _ 

12  FORMAT(/, ‘  FREQUENCY  ',F8.3,'  HZ*) 

C . TOE  NEXT  SECTION  OF  THE  PROGRAM  UP  TO  LINE  50  COMPUTES  THE  UPPER—. 

C  AND  LOWER  MODES  WICH  THE  PLATE  RESPONSE  IS  SUMMED  OVER.  N1  IS 

C.._  THE  LOWER  LIMIT  AND  N2  IS  THE  UPPER  LIMIT  _ 

C  INITIALIZE  AND  INCREMENT  N1 

_ Nlr-1 _ _ .  . .  _ 

18  Nl-Nl+1 

C  ._ CALCULATE  RESONANT  OMEGA  FOR  EN-Nl.  EM-N1  MODE  _ 

OMEGAMN0-(DA/(ROA»AaiA) )  •«  .  5*  ( ( (2 . 0«Nl+l .  0)  *PI/LX)**2. 0 

_ A*( (2.0*Nl+1.0)»PI/LY)**2. 0)  _  _ 

C  IF  RESONANT  OMEGA  IS  GREATER  THAN  OMEGA  GO  ON  TO  28  IF  NOT 

C _ INCREMENT_.N1  AND  REPEAT  _  _ _ 

IF  (OMEGAMNO.GT. OMEGA)  GOTO  28 

_ GOTO  .18  ... . . .  .  . 

C  INITIALIZE  AND  INCREMENT  N2 
OR  mo- 


38  N2-N2+1 

C-  CALCULATE  TOE  RESONANT  OMEGAS  FOR  EN=N2,  EMrQ  AND  EN-0,  EM— M2 
C  MODES  AND  TEST  IF  LARGER  THAN  OMEGA 

_ OMEGAMNL=(DA/(RQA*ACHA)  )•*  .5*.(  (  (2.Q?N2-tl-Q)!Pl/LX)*.*2 

t+(PI/LY)**2.0) 

_ OMEGAMN2?(DA/(ROA*ACHA)  )*?-.5!<(RI/LX)*?2.Q _ _ _ 

«+((2.0*N2+1.0)*PI/LY)**2.0) 

C _ IE  -LARGER  THAN  OMEGA  GO  ON  TO  48  IE  NOT  INCREMENT  N2  AND  REPEAT _ 

IF(  (OMEGAMN1.GT. OMEGA). OR.  (OMEGAMN2.GT. OMEGA))  GOTO  48 

_ GOTO  J8 _  _  _  _ 

48  CONTINUE 

C . ADO -THREE  MODES  TO  UPPER  LIMIT—  _ _ 

N2-N2+3 

_ N11?N] _ _  — . .  .  . . . 

C  SUBTRACT  THREE  MODES  FROM  LOWER  LIMIT,  ZERO  IS  LOWEST  N1 

_ N1-N1=3  _  .  .  .-.  .  _  .  .  _ 

IF(Nll.EQ.O)  Nl-0 

_ IE(Nll.EQ.l)  NlrO  .  _ _  _ 

IF(N11.EQ.2)  Nl-0 

_ WRITE. (6. 58)  N1.N2  .. .  ....  _  _ _ 

58  FORMAT  (215) 

C_ .THE  .PLATE  MODES  ARE  TAKEN  FROM  EM,  EN  =  N1  TO  N2.  _ 

DO  200  EM-N1.N2 

_ DO  200  EN-Nl. N2  ....  _ 

DO  300  I-l.N 

_ DO  300  J-l.N  ...  _ 

BT2(I,J)-BT(I,J) 

300  .  CONTINUE  _ 

C  OMEGAMNSQ  IS  THE  RESONANT  OMEGA* *2.0 

_ QMEGAMNSO-DSTAR? ( ( ( (2 . 0*EM+1 . 0)*P1/LX) **2.0f  ( (2.  O‘ENf-1.0) _ 

**PI/LY)»*2.0)**2.0)/(ROA*ACHA) 

-  .  DO  400  I-l.N  _  _ _ 

DO  400  J-l.N 

_ Xs(rLX/2. 0)  +  (LX*(  ItI)  )/N _  _ 

Y-C-LY/2.0)4-(LY»(J-l))/N 

_ B(  1 ,  J)?(COS(  (2. 0*£M+1 . 0)*P2*XO/LX)  )*(COS(  (2. 0*EN+1.  Q) _  . 

**PI*YO/LY))*(COS((2.0*EM+1.0)*PJ*X/LX))*(COS((2.0*ENU.O) 
**PI*Y/LY)  )/(OMEGAMNSQ-OMEGA**2 . 0) 


400 _ 


BT(I,J)-BT(I,J)+B(I,.J) 

CONTINUE 


200  CONTINUE 

_  WRITE  (26,301)  _ 

301  FORMAT  (*  PLATE  RESPONSE') 

_  DO  500  I-1,N  .  .  _ _ 

DO  500  J-l  ,N 

_  -  WA(l,J)-(4 . 0*F/MP) *BT{ I ,J) .  _ 

WAREAL(I,J)-WA(I,J) 

_ _ WA2(I , J )-(4 . 0*F/MP) *BT2( I, J)  .. _ 

WA2REAL ( I , J ) - WA2 ( I , J ) 

C_  WA  IS  TOE  COMPLEX  PLATE  RESPONSE  .  _  _ 

C  WAREAL  IS  TOE  REAL  PART  OF  THE  PLATE  RESPONSE 
C  WA2  IS  TOE  COMPLEX  PLATE  RESPONSE  OF  THE  PREVIOUS  MODE 
C  WA2REAL  IS  THE  REAL  PART  OF  THE  PREVIOUS  MODE 

500  CONTINUE  . . . _ 

C  THE  FOLLOWING  DO  LOOP  ADDS  C.\  THE  FINAL  POINT  TO  MAKE 

C  THE  MATRIX  SYMETRICAL  .  _  _ 

DO  550  I-1,N 

_  WAR£AL(N+l,I)-0.0  .  _  _ 

WAREAL( I ,N+1 )-Q . 0 

_  WA2REAL(N+1 . I )-Q. 0  _ 

WA2REAL(I,NU)-0.0 

_550  CONTINUE  _ 

C  OUTPUT  PLATE  RESPONSE 

DO  800  I-l.N+l  -  _ _ 

DO  800  J-l , N+l 

_  ..  HRITE  (26,401)  I-1.J-1,WA2REAL(I,J) ,1-1. J^l. WAREAL ( I, J) _ 

401  FORMAT  ('  WA2REAL( ' , 14 , 1 , ' , 14 , 1 )  -  ',£12.4,'  WAREAL ( ' , 

414^  •  ( ■  /  14  1 1  j  _  •  ,E12. 4)  _ _ _  _ 

C  CALCULATe'tHE  MAXIMUM.  MINIMUM,  AND  INCREMENT  FOR  TOE  PLOT 

_  VMAX-AMAXl ( ABS( WAREAL ( I , J) ) , ABS(VMAX) )  _ _ 

VMIN— VMAX 

_ VTNC-VMAX/2 . 0  _ _ 

800  CONTINUE 

C  PLOTTER  CREATES  A  3-D  PLOT  OF  TOE  PLATE  RESPONSE  _  _ _ _ 

CALL  PLOTTER  (VMAX , VMIN, VI NC, WAREAL, WORK1 ,WORK2 ,N, 

_ (OMEGA. ACHA. ETAA. LX, LY)  - - 

C1000  CONTINUE 

END  , . .  - - - 

C 

_  SUBROUTINE  PLOTTER  (VMAX,  VMIN,  VI NC,  WARE  AL,WQRKl,WORK2,N,_.. 

(OMEGA , ACHA , ETAA , LX , LY ) 

C_ PLOTTER  USES  TOE  ARL  HIDE/TEMPLATE  PLOTTING  PACKAGE  TO, CREATE _ 

C  3-D  PLOTS  OF  THE  PLATE  RESPONSE 

REAL  TYMIN, TYMIN, TXMAX.TYMAX.TXINC, TYING  _  ,  .. 

(OMEGA . ACHA , ETAA , PI , LX , LY 

_  PI-2. 0»ASIN( 1.0)  - . . . 

TXMIN-- LX/2 . 0 

TYMIN— LY/2.0  .....  _ 

TXMAX-LX/2 . 0 
,  ..  TYMAX-LY/2.0 
TXINC-LX/8.0 

TYINC-LY/0.0  .  _ 

WRITE  (6,501)  TXMIN,TXMAX,TXINC, TYMIN, TYMAX, 

(TYINC.N 

501  FORMAT  ( 3F12. 4 ,/, 3F12 . 4 , 14 ) 

CALL  USLPDF  . 

CALL  UPSET (  'OUTPUTFILE* ,7.0) 

CALL  UASSGN( 7.0,  ’i’RA.PDFV) 

CALL  USTART 

_ _  CALL  UDIM£N(  9. 0,6. 5)  .  _ _ 

CALL  USET( ’PERCENT’) 

CALL  UFONT(’CROM’) 

CALL  USET( 'LARGE') 

CALL  UPRINI’(61.  ,83.  , ’FREO-\’) 


U  H 


CALL  UPRI/ri  (OMEGA/ (2.0*  PI) ,  'REAL'  ) 

_ -CALL.  UPRINT(  77-  .  83 .  , ! OMEGA- \ '  ) _ _ 

CALL  UPRNT1  (OMEGA, 'REAL') 

_ CALL.UPRItfr(61.,73. ,  *H  -\!) _ _ 

CALL  UPRtfTl  ( ACHA ,  1  REAL  ‘ ) 

_ CALL  UPRIOT(77.,73., 'ETA-  V) _ . 

CALL  UPRWTI ( BTAA REAL 1 ) 

_ CALL-ARL  HLD£_FONT _ . 

CALL  UTOtfr('CROM') 

_ CALL  US£T(' LARGE'.) _ _ _ _ 

CALL  ARL  HIDE  SCALE(N+1,N*1,N+1,4000,N+1,VMIN,VMAX, 

_ A.'CUTOFF',cDtQFE7CL5,0.£,  1.25,0.5,0. 5, 0.0,45.0,45.0,  _ 

4'FDfiAW’,  'SDRAW') 

_ CALL-ARL  HIDE  <  HARE AL , W3RK1 , W0RK2 ) _ _ 

CALL  ARLllIDE  FAXIS( 'VIEW' , 1  REAL' ,TXMIN,TXItAX,TXINC, 

_ A' (F8. 2) I, 'X  .  METERS' ,1,50.0,50.0)  . 

CALL  ARL  HIDE  SAXIS(  'VIEW  ,  ‘REAL'  .TXMIN.TYMAX.'mNC, 

_ AM F8.2)',:5  METERS' ,1.50.0,50.0) 

CALL  ARL_HIDE_VAXIS( 'VIEW' , 'REAL' ,VMIN,VMAX,VINC, 

_  A' (E12. 4) ! .'AMPLITUDE' ,1.50.0,50.0) 

CALL  UEND 

_ RETURN _  __  _ _ 


END 
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_ PROGRAM  PR 

C_  DETERMINES. .THE  PLATE  RESPONSE  OF.  PLATE -B  USING  LMSL 
C  FFT  SUBROUTINE 

C _ _ _ _  _  _____ 

C  WRITTEN  BY  MICHAEL  F.  SHAW 

C _ _ _  _  _  _ 

INTEGER  OMEGAl,  O^GAB,  OMEGAT,  INC,  EM,  EN,  N,  I,  j, 

- &NROOEF , NCOOEF , LDCOEF , LDA , NLO , NUP , Nil  _ _ 

PARAMETER  (N-I28) 

- COMPLEX -DSTAR,OMEGAMNSO,COEF(N,N),A(N,N).Al(N,N),  _ 

4AT(N,N) ,KPF,KC,W(N,N) ,W1(N,N) , AT1(N,N) ,WC(N,N) , 

_ 4TB1P . . . . . 

REAL  WFF1(4«N+15),WFF2(4*N+15),CPY(N,N) 

-  REAL  VMAX,VMIN,VINC,WREAL(N+1.N+1),WORK1(4000),WORK2( 100) 

REAL  DA,DB,ETAA,ETAB,LX,LY,ROA,ROB, ACHA,ACHB,KX(N) ,KY(N) , 

- 4KX1 , KYI , KXM , KYN , F , RO , MP , XO , YO , IMN , 2 , PRA , GAMMA ,MU , NU , C , E , 

4HKX , HKY , PI , K , KP ,  OMEGA , OMEGAl , OMEGA2 , OMEGAMN0 , 

- AOMEGAMNl , OMEGAMN2 , W1REAL ( N+l.  N+l 3 , TEMP2 , TEMPI _ 

INTRINSIC  CMPLX 

- .EXTERNAL  F2T2B,FFTCI _  _ _ _  _ 

C  THE  NEXT  TWO  LINES  ALLOW  N  TO  BE  LARGER  THAN  64  POINTS 

- COMMON  /VORKSP/  RWKSP _ _ _ _ _ 

REAL  RWKSP (8084) 

c  VARIABLE  LIST 

C _  _ _ _ _  _ _ 

C  DA  -  BENDING  RIGIDITY  OF  TOE  FINITE  PLATE 

C  —  DB .-BENDING  RIGIDITY  OF  TOE  INFINITE  PLATE _ _ 

C  ETAA  -  DAMPING  COEFICIENT  FOR  TOE  FINITE  PLATE 

C — ETAB  -  DAMPING  COEFICIENT  FOR  TOE.  INFINITE  PLATE _ 

C  LX  4  LY  -  DIMENSIONS  OF  TOE  FINITE  PLATE 

C  —  ROA  DENSITY  OF  THE  FINITE  PLATE— _ 

C  ROB  -  DENSITY  OF  THE  INFINITE  PLATE 

C - ACHA_=_ THICKNESS  OF  TOE  FINITE  PLATE _ 

C  ACHB  -  THICKNESS  OF  TOE  INFINITE  PLATE 

C - RO  -  DENSITY  OF  FLUID  _ _  ...  _ _ 

C  XO  *  YO  -  TOE  POINT  WHERE  THE  INPUT  FORCE  F  IS  APPLIED 

C - Z  -  DISTANCE  BETWEEN  TOE  PLATES  _ _ _ _ 

C  PRA  -  PRANDEL  NUMBER  FOR  FLUID 

C — GAMMA  -  SPECIFIC  HEAT  RATIO  FOR  FLUID _  ...  _  . . . . 

C  MU  -  VISCOSITY  OF  FLUID 

C — NU  -  POISON'S  RATIO _  _ _  _  _  _  _ 

C  C  -  SPEED  OF  SOUND  IN  FLUID 

C - EA  .-.SHEER  MODULUS  OF  TOE  FINITE  PLATE  _ _ 

C  EB  -  SHEER  MODULUS  OF  TOE  INFINITE  PLATE 

C„_  101_r_ UPPER  WAVENUMBER  LIMIT  IN  FINITE  PLATE  IN  X  DIR  . 

C  KYI  -  UPPER  WAVENUMBER  LIMIT  IN  FINITE  PLATE  IN  Y  DIR 
C _ 

PARAMETER  (F-I.O) 

C -  INPUT  FORCE  TO  THE  FINITE  PLATE  IS  ASSUMED  TO  BE  A  UNIT  FORCE  (1) 

OPEN  (UNIT-25,  FILE- 'PR. DAT' ,  STATUS- ' OLD ' ) 

—  OPEN  (UNIT-26 ,  FILE- 'PR. OUT' .  STATUS-' NEW' ) 

C  ANOTHER  LINE  SO  N  CAN  BE  LARGER  THAN  64 

_  .  CALL. IWKIN( 8884) 

C  READ  INPUT  DATA  FROM  PR. DAT 

_ READ. (25,15)  LX.LY, XO, YO.AQ1A, ACHB 

IS  PORMAT(4F8. 3.2F8.6) 

- READ.  (25,25)  ROA, ROB, RO. ETAA, ETAB, GAMMA, PRA. _  _ 

25  F0RMAT(2F8.2,F8.3,2F8.7,2F8.6) 

READ. (25. 35)  MU.C.EA.EB.KXl.KYl  _ _ _  _ 

FORMATf  FB . 7 , F8 . 3 , 2E8 . 3 , 2F8 . 2 ) 

READ. (25, 45)  NU, OMEGAl, OMEGA2, INC. Z, OMEGA  _ 

FORMATf  F8 .4,F8.1,F16.1,I8,F8.3,F12.3) 

C  .COMPUTE  TOE  VALUE  OF  PI  _ _  ... 

PI-2. 0*ASIN( 1.0) 

C  COMPUTE  TOE  BENDING  RIGIDITIES 


nnonooo 
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DA-(EA*AaiA**3.0)/(]2.0*(  1 ,0-NU**2.0)) 

_ _  DB“(EB*ACHB**3 . 0)/( 12 . 0* ( 1. 0-NU**2. 0) ) 

DSTAR-CMPLX ( DA , DA* ( ETAA ) ) 

..  .  _ MP-ROA*ACHA*LX»LY  ..  .  .... 

C  THE  SAMPLE  SPACING  IN  WAVENUMBER  SPACE  IS  IIKX  AND  HKY 

_ HKX-KX1/N 

HKY-KYl/N 

C  PARAMETERS  FOR.  FFT2B  INVERSE  FOURIER  TRANSFORM  SUBROUTINE 
NRCOEF-N 

_...  _  NCCOEF-N. _ 

LDCOEF-N 

LDA-N _ 

WRITE(26,11)  LX.LY 

II  FORMAT ( ‘  PARAMETERS  OF  FINITE  PLATE' , /,  *  DIMENSIONS  OF  PLATE' 

(LX)  ’ ,F8. 3 , ’  (LY)  ' , FB . 3 , '  METERS') 

WRITE (26 >21)  XO, YO 

21  FORMAT (  1  DRIVING  POINT  XO  * ,F8.3,*  YO  ’,F8.3,' 

&  METERS')  ..  _ 

WRITE(26,31)  F 

31_ . _  FORMAT (  !  DRIVING  FORCE  \F8.3,'  N') 

WRITE(26,41)ACHA 

.41 _ FORMAT ( !  PLATE  THICKNESS  '  ,F8.6,  '  M*). 

WRITE (26, 51)  ROA 

,  51. ..  FORMAT!  1  PLATE  DENSITY  ',F8.3,'  KG/CU-METERS ' ) 

WRITE(26,61)  EA 

61 _ FORMAT ( '  MODULUS  OF  ELASTICITY  '  ,E8.3,  *  N/SQ-METERS * ) 

WRITE(26,71)  DA 

_71 _ FORMAT( !  BENDING  RIGIDITY  *,E12.4)  _ 

WRITE! 26, 81)  ETAA 

81 _ FORMAT!'  DAMPING  COEFFICIENT*.  F8. 4  >  -  ...... 

WRITE(26,91) 

_91 _ FORMAT!/,*  PARAMETERS  OF  INFINITE  PLATE’./) 

WRITE (26, 101)  ACHB 

JL01 . .  FORMAT! 1  PLATE  THICKNESS  *,F8.6,'  M' >  - 

WRITE (26, 111)  ROB 

III  _ FORMAT!  *  PLATE  DENSITY  \F8.3,'  KG/CU-METERS')  ...  - 

WRITE(26,121)  EB 

-121 _  FORMAT! '.  MODULUS  OF  ELASTICITY  .’,E8.3,'  N/SQ-METERS  *  )  ...  . 

WRITE(26,131)  DB 

131 _ FORMAT!  1  BENDING  RIGIDITY  *,E12.4)  _  . . - 

WRITE(26 ,141)  ETAB 

141 _ FORMAT!!  DAMPING  COEFICIENT*  ,F8. 4)  - 

WRITE(26,151)  RO 

J3L _ FORMAT!/, '.  PARAMETERS  FOR  FLUID',/,*  DENSITY  '.F8.3,  - 


4*  KG/CU-METER  * ) 


..WRITE (26, 161)  MU 

FORMAT! *  VISCOSITY  ' ,F8. 

- ; 

16  i 

7,'  KG/M-S*) 

WRITE(26,181)  C  .  . 
FORMAT! '  SPEED  OF  SOUND 

. 

- - 

i8i 

' , F8 . 3 , '  M/S*) 

WRITE(26,191)  PRA 

. .  i 

191 

FORMAT!*  PRANDTL  NUMBER 
WRITE! 26, 201)  GAMMA 

* ,F8.3) 

201 

FORMAT!*  SPECIFIC  HEAT  RATIO  \F8.3) 

WRITE(26,211)  Z 

_  j 

211 

FORMAT!*  DISTANCE  BETWEEN 

TLATES  ' ,F8.3, '  M*) 

i 

C  INITIALISATION  ROUTINES  FOR  THE  FFT 
CALL  FFTCI(N.WFFl) 

_ CALL  FFTCI(N,WFF2)..  .  .  - -  -  ....  -  - 

THE  FREQUENCY  RANGE  IS  TAKEN  FROM  OMEGAB  TO  OMEGAT  IN  INCREMENTS 
_ CORRESPONDING  TO  INC  .  .  .  -  -  -  - 

_ OMEGAB"10*IFIX(ALOG10(OMEGA1) )  -  - 

OMEGA2-INC*10» I FIX ( ALOG10 (OMEGA2 ) ) 

DO  1000  OMEGAI-OMEGAD, OMEGAT  —  ■ 

OMEGA-2 .0'PI •  ( 10. 0»«  (OMEGAI/! INC»1 0. 0) ) ) 

DO  100  I-l.N 
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DO  100  J-1,N 

_ AT(I,J)-0.0 _  _  _ 

100  CONTINUE  ' 

C_  It  IS  THE  WAVENUMBER  IN  TOE  FLUID  — 

K^JMECA/C 

C  KC..IS  TOE  COMPLEX  WAVENUMBER _  _  _ 

KC- ( K )  *CMPLX  (1.0,(OMEGA*MU*(  (4. 6/3. 6  )•*■((  GAMMA-1. 0)/PRA) 

S)/(2.Q*RO*C**2.0) ) ) _ _  _  _ _  _ 

C  KP  IS  THE  PLATE  WAVENUNBER 

- KP-(ROB*ACHB*OM£GA**2.0/DB)**.25  _ _ 

C  KPF  IS  TOE  COMPLEX  WAVENUMBER**4 . 0 

- KPF-(ROB*ACHB*OMEGA**2.0/DB)*OlPLX(1.0,ETAB)  _  __  _  _ 

WRITE (26. 12)  OMEGA/(2.0*PI) 

12 - FORMATf/, 1  FFEOUENCX  ‘.F8.3.'  HZ') 

WRITE (26 ,22)  KC 

22 -  FORMAT ('  COMPLEX  WAVENUMBER  *.F8.3,'  +  i  ‘,F8.3,*  1/M') 

WRITE(26 , 32)  (KPF)**. 25 

-  -32—  FORMAT  {  '  COMP  LEX  _P  LATE  WAVENUMBER  ',E12.4,'  +  i  ' 

&.E12.4, '  1/M') 

C  _ THE  NEXT  SECTION  (UNTIL  .LINE. 58)  COMPUTES  TOE  UPPER  AND  LOWER _ 

C  MODES  WICH  TOE  EQUATION  IS  SUMMED  OVER.  NUP  IS  TOE  UPPER  LIMIT 

C_AND  NLO.IS.THE  LOWER.  LIMIT _ _ . _ _ 

C  INITIALIZE  AND  INCREMENT  NLO 

_ NLO-rl  _ _  _ _ _ _ 

18  NLO-NUHl 

C — CALCULATE  .RESONANT  OMEGA  FOR  TOE.  EN-NLO.  EM-NLO  MODE  _ _ 

OMEGAMN0-(DA/( ROA* ACHA)  )  ••  .  5*  ( ( ( 2 . 0*NIXH1 . 0)  *PI/LX)  «*2 . 0 

- 44(  (2.0*NIXHl.0)*PI/LY)**2.0) _ : _ 

C  IF  TOIS  RESONANT  OMEGA  IS  LARGER  THAN  OMEGA  GO  ON  TO  28  IF  NOT 

C _ INCREMENT  NLO  AND  REPEAT  _ _ _ _ _ _ 

IF  (OMEGAMNO.GT. OMEGA)  GOTO  28 

_ GOTO  .18  _ _ _ _ 

C  INITIALIZE  AND  INCREMENT  NUP 

.28 _ NUPrrl _  _ _ _ _ 

38  NUP-NUP+1 

C-  CALCULATE  TOE  RESONANT  OMEGAS  FOR  EN-NUP,  EM-0  AND  EN-0,  EM-NUP 
C  MODES  AND  TEST  IF  LARGER  TORN  OMEGA 

- OMEGAMNl-(DA/(ROA*AOiA)  ) ••  .5* ( ( (2 . 0*NUP+1. 0)  *PI/LX)  **2 .0 _ 

4+(PI/LX)**2 . 0) 

- OMEGAMN2-(DA/(ROA*ACHA)  )•• .  5*(  (PI/LX)**2.0  _ _ _ 

t* ( ( 2 . 0*NUP*1 . 0 ) *PI/LT ) • • 2 . 0 ) 

C.-XF. LARGER  THAN  OMEGAGO  ON  TO  48,  IF  NOT  INCREMENT  NUP  AND  REPEAT _ 

IF(  (OMEGAMN1.GT. OMEGA)  .OR.  (OMEGAMN2.GT. OMEGA))  GOTO  48 

_ GOTO  38  _ _ _ _ _ 

48  CONTINUE 

[C_  ADD.  TOREE  MODES  ON  TO  UPPER  LIMIT  ..  ......  _  . 

NUP-NUP+3 

.  ...  .  -N11-NL0  _  _  ...  .....  _ 

C  SUBTRACT  THREE  MODES  FROM  LOWER  LIMIT,  ZERO  IS  THE  LOWEST 

_ NLO- NLO- 3  _  _ 

IF(Nll.EO.O)  NLO-O 

_  IF(Nll.EO.l)  NLO-O  .  .  . 

IF(N11.E0.2)  NLO-O 

_  WRITE  (6,58)  NLO, NUP  . 

58  FORMAT  (215) 

C_TOE. PLATE  MOOES  ARE  TAKEN  FROM  EM  AND  EN  -  NLO  TO  NUP  _ 

DO  200  EM-NLO, NUP 

_ DO.  200  EN-NLO,  NUP _ _ 

DO  300  I— 1,N 

_ DO -300  J-l.N _ _ _ _ 

C  ATI  IS  TOE  PREVIOUS  MODE  SUM 

_ AT1(I,J)-AT(I,J) _ _  _ _ _ 

300  CONTINUE 

C  CMEGAMNSO  IS  TOE  RESONANT  OMEGA *2  ...  .  _ _ 

OMEGAMNSO-DSTAR*  (  (  (  ( 2 . 0*B1*  l .  0)  *PI/I.X) **2 . 0* ( (2 . 0*EN+1 . 0) 
••PI/LX)  **2.0)**2.0)/(  ROA*  ACHA) 


UUUUOUUUUU  I  '  I  I  OO 


KXM-((2.0»EM*1.0)«PI/I.X) 

KYN«?((2.0«ENU.0)*ri/LY) 

C  COEFICIENT  CALCULATES  THE  INPUT  ARRAY  (COEF)  TO  THE  FFT 
.  CALL  COEFICIENT  (N,HKX,HKY.KXl,m.PI  .KXM.KYN. 
«K,KC,KPF,Z,KX,KY,LX,LY,EM,nN,COEF) 

C  -  F2T2B  IS  AN  IMSL  SUBROUTINE  MUCH  CALCULATES  THE  2-D  FFT  OF 
C  A  SET  OF  FOURIER  COEFFICIENTS. 

C  COEF( I ,  J)  IS  HIE  INPUT  ARRAY 
C  A( I , J )  IS  THE  OUTPUT  ARRAY 

.  _  ....  CALL  F2T28  ( NR  COEF ,  NCCOEF ,  A ,  LDA ,  COEF ,  LDCOEF ,  WFF1 .  WFF2 , 

tWC,CPY) 

.  , . .  DO  400  I-1,N 

DO  400  J-1,N 

..  A1(I, J)-( 1/(2. 0*PI)**2. 0)  *I1KX*IIKY*  ( (-1)  •*!  )  *  ( (-1) **J) 

4*A(I,J)'(COS((2.0»EM+1.0)*P!*XO/LX)* 

.  6 COS ( ( 2 . 0  *  EN+ 1 . 0 ) • P I *  YO/L Y ) / { OMEG AMNSQ-OMEGA*  *  2.0)) 

AT(I , J)-AT( I , J)+A1( I , J) 

400  CONTINUE 

200  CONTINUE 

_ WRITE  (26,301). 

301  FORMAT  ('  PLATE  RESPONSE') 

_ DO  500  Irl.N  .  ^ 

DO  500  J-1,N^5 

_  W(  I  ,J)“(  *F*RO*OM£GA**2. 0*  AT(  I ,  J) )/(MP*DB) 

Wl(I,J)-( (-4.0),F»RO«OMEGA»«2.0»AT1(I.J) )/(MP«DB) 

. .  WREALt I , J ) -W(I 7J 7s  & 

W1REAL(I,J)-W1(I,J) 

C_  WREAL(I.J)  IS  THE  REAL  PART  OF  THE  PLATE  RESPONSE 
C  W(I,J)  IS  THE  COMPLEX  PLATE  RESPONSE 

C W1(I,  J)_  ......IS  THE  COMPLEX  PLATE  RESPONSE  OF  THE  PREVIOUS  MODE 

500  CONTINUE 

_THE_FOLLOWING. TWO  DO-LOOPS  CONVERT  THE  OUTPUT  PLATE  RESPONSE 
IN  THE  X-Y  PLANE  SO  THAT  THE  RESPONSE  SHOWS  THE  CENTER  OF  THE 
_ PLATE  IN  THE  CENTER  OF  THE  ARRAY 
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_  .  DO  600  I-1,N 

DO  6P0  J-l,N/2 
TEMPl-WREAL(I.J) 

TEMP2-W1REAL(I,J) 

WREAL( I , J ) -WREAL ( I , N/2 d ) 

W1REAL ( 1 , J ) -W1 RE AL ( I , N/2  + J ) 

WREAL( I , N/2+ J ) -TEMT1 
W1REAL  ( I ,  N/2-t  J )  -TEMP2 
500  CONTINUE 

DO  700  J-l.N 
DO  700  I-l.N/2 
TEMPI “KREAL(  I  ,.I) 

TEMP2-W1REAL(I.J) 

WREAL( I , J ) -WREAL ( N/2 i I , J ) 

_ _ WlREAL(I.J)-WlREAL(N/2+I.J) 

WREAL( N/2+I , J )-TEMPl 

_ WlREAL(N/2+I,J)-TEMP2 

(00  CONTINUE 

_T»IE  FOLLOWING  DO  LOOP  ADOS  ON  THE  FINAL  POINT  TO  MAKE 
THE  MATRIX  SYMETTUCAL 

_  .  DO  550  I-1,N»1 

WREALt I ,N+1 ) -WREAL ( 1,1) 

WREAL ( N* 1 . I ) -WREAL ( 1 , I ) 

550  CONTINUE 

C  OUTPUT  PLATE  RESPONSE _ 


DO  800  I«1,N 

_ DO  800  J-l.N _ _  _ _ 

C  WRITE  (26,401)  I-l,J-l,WlREAL(i,j),i-l,J-l'wREAL(f,J) 

C  401 - FORMAT .(!.  Wl( !,  14 I4.J  ).  -  1  ,E12.4,2X _ 

C  414,  V, 14,')  -  '  ,£12.4) 

C ..CALCULATE. THE  MAXIMUM,  MINIMUM.  AND  INCREMENT  FOR  TOE  PLOT _ 

VMAX-AMAX1 ( ABS ( WREAL ( I , J ) ) , ABS ( VMAX ) ) 

- WilNr-VMAX _  _  _  _ 

VlNC-VMAX/2.0 

800. _ CONTINUE  .  _  .  _  _  _  _ _ 

C  PLOTTER  CREATES  A  3-D  GRAPH  OF  THE  PLATE  RESPONSE 

- CALL-PLOTTER  ( VMAX, VMIN.VINC, WREAL, WORKl,WORK2,N,KXl, KYI,  _ 

4 2 , OMEGA , ACHB , CTAB , KP , K ) 

C  1000  CONTINUE  _.  _ 

END 

C  .....  _ 

SUBROUTINE  COEFICIENT  (N,HKX,HKY,KX1,KY1,PI,KXM,KYN, 
4K.KC,KPF.Z,KX,KY,LX,LY,EM,EN,COEF)  .  - 

C  COEFICIENT  CALCULATES  THE  FOURIER  COEFICIENTS  FOR  USE  IN  THE 

C  FFT2B  SUBROUTINE  . .  ...  .....  .  _ _  _ 

REAL  KX(N),KY{N),IMN,HKX,HKY,PI,KXl(m, 

- 4KXM.KYN,Z.LX,LY.X . . . . . .  . 

INTEGER  EM,EN,N 

_ COMPLEX  TEMP,KC,KPF,COEF(N,N) . . .  .  . 

DO  10  I-l.N 

_ DO  _10- J“1.N  _ _ _  _ _  _ _ 

C  THE  FFT  USES  THE  INTERVAL  FROM  -KX1/2  0  TO  KX1/2.0 

- KX(I)-((I-l)«HKX)-(KXl/2.0).  _ _ _ _ _ 

KY(J)-((J-l)«HKY)-(KYl/2.0) 

C  TOE.  FOLLOWING  IF  STATEMENTS  DETERMINE  THE  VALUE  FOR  IMN  — _ 

IF  <<XXM**2.0  .EQ.  KX(I)**2.0) .AND. 

_ 4(KYN*«2-0..EQ.  KY(J)**2.0) ) _ _ 

4IMN-LX-LY/4.0 

_ IE.  (  (KXM**2.0  . EQ.  . KX(I)**2.0)  .AND.  .  _ _ 

4(KYN**2.0  .NE.  KY(J)**2.0)) 

_ 4IMN-LX«KYN«COS(KY(J)*LY/2.0)R(-1.0)**EN/  _  _ _ 

4(KYN**2.0-KY( J) **2. 0) 

_ IE_((KXM««2.0  .NE.,KX(I)*«2.0).AND.  .  ...  _  _ 

4(KYN*«2.0  .BO.  KY(J)**2.0) ) 

_ 4  IMN~LY*KXM*OOS  ( KX  ( I )  *LX/2 . 0  )  •  (“1 . 0 )  *  *EM/ _ 

4(KXM**2.0-KX(I)*«2.0) 

_ IF_((KXM**2.0  .NE.-KX(I)**2.0).AND.  . . .  .  _ _ _ _ 

4(KYN**2.0  .NE.  KY(J)»«2.0)) 

_ 4lMN-4-0»KXM*KYN*COS(KX(I)*LX/2.0)«  ......  _ _ 

4COS(KY(J)*LY/2.0)«(-1.0)»*EM*(-10)»*EN/ 

....  4{(KXM««2.0-KX{I)»«2.0)*(KYN«*2.0-  ....  ...  .  . . . 

4KY( J)**2. 0)  ) 

C  TOE  FOLLOWING  IF-STATEMENT  TESTS  IF  TOE  EXPRESSION  IS  ......  _ 

C  NEGATIVE  FOR  USE  IN  TOE  COMPLEX  EXPRESSION  FOR  COEF(I,J) 

_ IF  (K*»2.0-KX(I)**2.0-KY(J)««2.0)  55.  65,  65  ....  ..  _ 

55  TEMP-01PLX({{ABS(K»«2.0-KX(I)*«2.0- 

.  4KY(J)»*2.0))*«.5)«(-Z),0.0) 

GOTO  75 

65  TEMP-CMPLX(0.0,((K«*2.0-KX(I)*'2.0- 

4KY(J)**2.0)««0.5)»Z) 

75-  ...  CONTINUE  _ _  ... 

C  CALCULATE  OOEF 

_ COEF{I.J)"CEXP(TQ!P)*IMN/_ _  .  .  ..  _ _  _ 

4( (KC**2.0-KX( I) ••2.0-KY(J),*2.0) •• .  5* 

_ 4((KX(I)*«2.0+KY(J)*«2.0)»»2.0-KPF))  _ _ _ _ 

10  CONTINUE 

_ RETURN _ _ .  _ _ 


SUBROUTINE  PJ/TTTEF  (VMAX,VMIN,VJNC,WREAL,NORKl ,WT>RK2,N, 
4KX1 , KYI , Z .OMEGA , AQID, CTAB ,  KT , K) 


C  PIOTTER  USES  THE  APLJUDC/TCIPLATE  PLOTTING  PACKAGE  TO  CREATE 
C_3-D  PLOTS  OF . THE  PLATE  RESPONSE 

REAL  KX1, KYI, TXMIN, TYMIN, TXMAX,TYMAX,TXINC,TYINC 

- 4QMEGA,  AQ1B,ETAB,KP,K,  Z,PI  .  .. 

PI-2.0*ASIN{1.0) 

_ TXMINr-N«PI/KXl  ....  . 

TYMIN— N«PI/KY1 
_  ....  _  TXMAX-N«PI/KX1  ... 

TYMAX-N»PI/KYl 

— . .  TXINC-N*PI/{ 4 . 0*KX1 ) 

TYINC-N*PI/{ 4 . 0*KY1 > 

...  - WRITE  (6,501)  JCXl, TXMIN, TXMAX.TXINC, KYI, TYMIN. TYMAX, 

4TYINC.N 

501  FORMAT  (4F12.4,/.4F12.4,I4) 

CALL  USLPDF 

CALL  UPSET(  'OUTPUTFILE ',7.0) 

CALL  UASSGN(7 . 0,  'PR. PDFS') 

_  CALL  USTART.  _ 

CALL  UDIMEN(  9.0,6. 5) 

_ CALL  USETC PERCENT' ) 

CALL  UPONT('CROM') 

_ CALL  USET( 'LARGE' ) 

CALL  UPRIWT(61. ,93. , 'KX1-  \') 

_ CALL  UPRNT1(KX1, 'REAL' )  __ 

CALL  UPRINT(77. ,93. , *KYl-  \') 

_ CALL  UPRNT1(KY1, ‘REAL* ) 

CALL  UPRINT(61. ,88. , 'Z  -  \' ) 

_ CALL  UPRNTl ( Z , 'REAL' )  . 

CALL  UPRINT(77. ,88. , 'H  -  \' ) 

_ _ CALL: UPRNTl ( ACHB ,  'REAL'  )  . 

CALL  UPRIMTf  61 . ,  83 .  ,  '  FREO-\ '  ). 

_ CALL  UPRNTl(OMEGA/(2. 0*PI) , 'REAL* ) 

CALL  UPRINT( 77 . , 83 . , ’ OMEGA-\ * ) 

_ CALL  UPRNTl (OMEX3 A,  'REAL' )  ..  .  _ 

CALL  UPRINT(61. ,78. , ’K  -\’) 

_ CALL  UPRNTKK, ’.REAL')  .  _ 

CALL  UPRINT(77. ,78. , ’ KP  -  \’) 

_ CALL. UPRNI1(KE, ’REAL'  ) 

CALL  UPRINT(77. ,73. , 'ETA-  \’) 

_  CALL  UPRNTl (ETAB, 'REAL' ) 

CALL  ARL_HIDE  FONT 

_ CALL  UFONT( .'  CROM ' )  . 

CALL  USET( 'LARGE' ) 

_ CALL.ARL  HIDE  SCALE(N+1  ,N+1  ,N+1,4000,NU,VMIN,VMAX,  _  ... 

A 'CUTOFF' ,CUTOFF70. 5, 0.6, 1.25, 0.5, 0.5, 0.0, 45. 0,45.0, 

.....  4'FDRAW' ,  'SDRAW' ) _ 

CALL  ARL  HIDE(WREAL,WORKl,WORK2) 

_ CALL  ARL~HIDE  FAXIS(  'VIEW'  ,  'REAL'  ,TXMIN,TXMAX.TXINC, 

4' (F8.2) ' , 'X  METERS' ,1,50.0,50.0) 

....  ...CALL  ARL  HIDE  SAXISC  'VIEW1 ,  'REAL'  .TYMIN,  TYMAX.TYINC, 

4 ' { F8 . 2 ) ' , ' Y  METERS' ,1,50.0,50.0) 

.  CALL  ARL  HIDE  VAXISC  'VIEW'  ,  'REAL'  ,VMIN,VMAX, VINC, 

4' (E12.4) ' ,TAMPLlTUDE' ,1,50.0,50.0) 

CALL  UEND 
RETURN 

_ END- _  .  .  . 
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PROGRAM  PRTL5K 

.  DETERMINES  THE  TRANSMISSION  LOSS  BETWEEN  PLATE-A  AND  PLATE-B  .  _ 

USING  RMS  VIBRATION  LEVELS  OF  EACH  PLATE 

WRITTEN  BY  MICHAEL  F.  SHAW 

INTEGER  OMEGAI ,  OMEGAB,  OMEGAT,  INC,  EM,  EN,  I,  J, 

_  ANRCOEF ,  NCCOEF ,  LDCOEF ,  LDA ,  N ,  NIO ,  NUP ,  Nl  1  _ _ 

PARAMETER  (N-128) 

COMPLEX  DSTAR , OMEGAMNSQ , COEF { N , N ) ,A(N,N),B1(N,N). 
4BT(N,N),KPF,KC,WA(N,N),WB(N,N),AT{N,N),B(N,N), 

AWC(N.N)  ...  ..  - . 

REAL  WFF1(4«N+15) ,WFF2  ( 4«N<-15) ,CPY(N,N) 

REAL  WBREAL(N.N) ,WAREAL(N,N) ,WORKl(4000) ,WORK2(100) , 
4AMPA,AMPB,TAU,TL(200) , COUNT, OMEGAI  ,OMEGA2,  Xl  (200) , OMEGART 
REAL  DA ,  DB ,  ETAA ,  ETAB ,  LX ,  LX ,  ROA ,  ROU ,  ACi  I A ,  ACHB , KX ( N ) , K Y ( N ) , 
4KXl,KYl,KXM,KYN,F,RO,MP,XO,YO,IMN,2,PRA,GAMMA,MU,NU,C,E, 

_  41  1KX . HKY , P I ,  K ,  KP , OMEGA , OMEG AMN 0 , OMEG AMN 1 , OMEGAMN2  .  .  ... 

INTRINSIC  CMPLX 

_  EXTERNAL  F2T2B.FFTCI  - - -  . 

C  THE  NEXT  TWO  LINES  ALLOW  N  TO  BE  LARGER  THAN  64  POINTS 

_  COMMON  /WORKSP/  RWKSP  _ _ _ 

REAL  RWKSP (8884) 

C  ... 

C  VARIABLE  LIST 

C  DA  -  BENDING~RIGIDITY  OF  TOE  FINITE  PLATE 

C  DB  -BENDING  RIGIDITY  OF  THE  INFINITE  PLATE  _ _  _ 

C  ETAA  -  DAMPING  COEFICIENT  FOR  THE  FINITE  PLATE 

C  ... . ETAB  -  DAMPING  COEFICIENT  FOR  THE  INFINITE  PLATE  -  - 

C  LX  A  LY  -  DIMENSIONS  OF  TOE  FINITE  PLATE 

C  _  ROA  ~  DENSITY  OF  THE  FINITE  PLATE  -  _ ....  —  - 

C  ROB  -  DENSITY  OF  THE  INFINITE  PLATE 

C._  ACHA  -  THICKNESS  OF  THE  FINITE  PLATE  -  - 

C  ACHB  -  THICKNESS  OF  THE  INFINITE  PLATE 

C  _  RO  -  DENSITY  OF  FLUID  -  - - - 

C  XO  4  YO  -  TOE  POINT  WHERE  TOE  INPUT  FORCE  F  IS  APPLIED 

C  2  -  DISTANCE  BETWEEN  THE  PLATES  - 

C  PRA  -  PRANDEL  NUMBER  FOR  FLUID 

C  _  GAMMA  -  SPECIFIC  HEAT  RATIO  FOR  FLUID  -  —  — 

C  MU  -  VISCOSITY  OF  FLUID 

C_  NU  -  POISON'S  RATIO  --  - 

C  C  -  SPEED  OF  SOUND  IN  FLUID 

C_  EA  -  SHEER  MODULUS  OF  THE  FINITE  PLATE  ..  - 

C  EB  -  SHEER  MODULUS  OF  TOE  INFINITE  PLATE 
C  KXl  -  UPPER  WAVENUMBER  LIMIT  IN  FINITE  PLATE  IN  X  DIR  _ 

C  KYI  -  UPPER  WAVENUMBER  LIMIT  IN  FINITE  PIATE  IN  Y  DIR 

C_. 

PARAMETER  (F-1.0) 

C  INPUT  FORCE  TO  THE  FINITE  PLATE  IS  ASSUMED  TO  BE  A  UNIT  FORCE  (I) 
OPEN  (UNIT-25,  FILE-' PR5K. DAT* ,  STATUS- 'OLD' ) 

OPEN  (UNIT-26,  FILE- ' PRTL5K. OUT' ,  STATUS- 'NEW' ) 

C  ANOTHER  LINE  SO  N  CAN  BE  LARGER  THAN  64 
[  _  CALL  IWKIN(8884) 

C  READ  INPUT  DATA  FROM  PR. DAT 

READ  (25,15)  LX ,  LY ,  XO ,  YO ,  AQIA ,  ACI  111 
15  PORMAT( 4F8. 3.2F8.6) 

...  .  READ  (25.25)  ROA. ROB, RO, ETAA, ETAU . GAMMA , PUA 

25  FORMAT ( 2F8 . 2 , F8 . 3 , 2F8 . 7 , 2F8 . 6 ) 

READ  (25.35)  MU.C.EA.ED.KXl.KYl 
35  FORMAT(F8.7,F8.3,2E8. 3.2F0.2) 

READ  (25,45)  NU. OMEGAI. OMEGA2. INC, 2, OMEGA  .  .. 

45  FORMAT( F8. 4 ,F8.1 ,F16. 1 , I0,F8. 3.F12. 3) 

C  COMPUTE  THE  VALUE  OF  PI 
PI-2. 0«ASIN( 1.0) 

C  COMPUTE  THE  BENDING  RIGIDITIES  _ 
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DA-{EA*ACHA«*3 . 0)/( 12 . 0« ( 1 . 0-NU««2 . 0) ) 

_ DB-(EB*ACHB,?3. Q)/(12 . 0*  ( 1.0~NU**2. 0) ) _ _ 

DSTAR-CMPLX ( DA , DA*ETAA) 

_  MP-R0A*ACJ1A*LX«LY _  _  _ _ _ 

WRITE (26, 11 )  LX, lX 

11- .  _  FORMAT ( *  PARAMETERS  OF  FINITE  PLATE*,/,!  DIMENSIONS  OP  PLATE' 
4,*  (LX)  1 ,F8.3, 1  (LY)  ‘ , F8 . 3 , '  METERS') 

_ WRITE(26,21)  XO.YO  .  _  .  _ 

21  FORMAT ( *  DRIVING  POINT  XO  ',F8.3,'  YO  *,F8.3,' 

_ 4  METERS')  ..  _  _  _  _ 

WRITE(26,31)  F 

-31 _ FORMAT { DRIVING  FORCE  .  ' ,F8. 3,1_.  N! ) _  .  _ _ 

WRITE(26 ,41 ) ACHA 

41 _ FORMAT  ( !  PLATE  THICKNESS  ', F8.fi,!  M‘ )  _  _ _ _ 

WRITE (26, 51)  ROA 

-51  .  FORMAT ( '  PLATE  DENSITY  '  ,F8. 3,1_  KG/CtJ-METERS' )  _ 

WRITE(26,61)  EA 

61 _ FORMAT  ( !-MODULUS  OF  £LASTICITY.J-,E8.2,  '  N/SQ-METERS '. ) _ 

WRITE(26,71)  DA 

_21 _ FORMAT(  L  BENDING  RIGIDITY  !,E12.4).._. _  _ 

WRITE (26, 81)  ETAA 

-81 _ PORMATC  DAMPING  COEFFICIENT'  ,F8.4)_  _  _ 

WRITE(26,91) 

91 _ FORMAT(/-,  !  PARAMETERS  OF  INFINITE  PLATE'  ,/)  _ 

WRITE (26 ,101)  ACHB 

101  FORMAT (  L. PLATE  THICKNESS  -  1  ,F8.fi^l_  M'  ) _  .  _ 

WRITE (26, 111)  ROB 

111 _ F0RMAT( ‘.-PLATE  DENSITY _ L  ,E8.3,1_  KG/CU-METERS ' ) _ 

WRITE (26, 121)  EB 

-121 _ FORMAT ( ! -MODULUS  OF.  ELASTICITY !^E8.3,  !.  N/SQ-METERS ' ) 

WRITE(26,131)  DB 

131 _ FORMAT (  L  BENDING. RIGIDITY — '  ,E12.4). _  _ 

WRITE(26 ,141)  ETAB 

-141 _ FORMAT ( — . DAMPING  (DEFICIENT  ,F8. 4) _ _ _ 

WRITE (26, 151)  RO 

-151 _ FORMAT  ( / , '  PARAMETERS  FOR  FLUID!,/.,!  DENSITY  \F8.3, _ 

4'  KG/CU-METER') 

_ HRITE(26,161)  MU . . . . .  .  _ 

161  PORMATC  VISCOSITY  ',F8.7,'  KG/M-S') 

_ WRITE (26, 181)  C  _  -  . . .  . 

181  PORMATC  SPEED  OF  SOUND  ',F8.3,'  M/S') 

_ WRITE (26, 191)  PRA—  _  . _ 

191  PORMATC  PRANDTL  NUMBER  ',F8.3) 

_ WRITE  (26, 201)  GAMMA _  _ _ _ _ 

201  PORMATC  SPECIFIC  HEAT  RATIO  ’,F8.3) 

_ WRITE(26,211)  Z  _  ..  _ 

211  PORMATC  DISTANCE  BETWEEN  PLATES  \FB.3,'  M‘) 

C  TOE  FREOUENCY  RANGe"  IS~ TAKI^lROM  1258  TO '  5000  HZ  ™ 

_ DQ  1000  QMEGAI-31,37  .  _ 

OMEGA-2 . 0»PI« ( 10 . 0»« (OMEGAI/10. 0) ) 

C  _  KX1  AND  KYI  ARE . .  . 

C  THE  LIMITS  OF  THE  FFT.  THE  FFT  IS  TAKEN  FROM  -KXl/2.0  TO 

C  KK1/2.0  AND  FROM  rKYl/2.0  TO  KYl/2.0  .  _ 

KX 1-4 02. 2 

_ KYl-4  02. 2 _  _  _ _ 

C  THE  SAMPLE  SPACING  IS  HKX  AND  HKY 

_ HKX-KX1/N _  _  -  . .  . 

HKY-KY1/N 

C_  PARAMETERS  FOR.FFT2B  INVERSE  FOURIER  TRANSFORM  SUBROUTINE  _ 

NRCOEF-N 

_ NCCDEE-N . .  . .  . . 

LDOOEF-N 

_ _ LDA-N _ _ _  _  _  _ 

C  INITIALIZATION  ROUTINES  FOR  THE  FFT 
CALL  FFTCI (N, WFF1) 


nnnnnnnnnn 
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CALL  FFTC1 (N.WFF2) 

C_  lt.IS  THE  WAVENUMBER  IN  THE  FLUID  - - 

K-OMEGA/C 

C  KC.IS  THE  COMPLEX  WAVENUMBER  -  - 

KC- ( K ) *CMPLX (1.0,  ( OMEGA*MU«(( 4. 0/3. 0)+( (GAMMA-1. 0)/PRA) 

i)/(2. 0*RO*C**2. 0) ) )  -  -  -  - 

C  KP  IS  THE  PLATE  WAVENUMBER 

KP-(ROB*ACHB*QM£GA**2 . 0/DB) ••  .25  -  - - - 

C  KPF  IS  THE  COMP1.EX  PLATE  WAVENUMBER»«4 . 0 

.  KPF“(ROB*AQIB*OMEGA**2 . 0/DU)  *CMPLX(  1 . 0 ,ETAB)  . . 

DO  100  I-l.N 

_  DO  100  J-l.N  ..  ..  - 

BT( I , J)-0. 0 

.... _ AT(I,J)-0.0  .  -  -  - .  ■  - — 

100  CONTINUE 

C  QMEGAMNSUB  CALCULATES  THE  UPPER  AND  LOWER  MODES  NUP,  NLQ  .  ...  - 

CALL  OMEGAMNSUB ( NLO , NUP , OMEGA , RO A , AC1 1 A , DA , LX , LY ) 

C._.  THE  PLATE  MODES  ARE  TAKEN  FROM  EM  AND  EN  -  NLO  TO  NUP  -  ... - 

DO  200  EM-NLO,NUP 

_ ...._.  DO  200  EN-NLO.NUP  -  - 

C  OMEGANMSQ  IS  THE  RESONANT  OMEGA* *2 

_  QMEGAMNSQ-DSTAR* ( ( ( (2 . 0*EM+1. 0) *PI/LX) **2 . 0+ ( (2.0*EN+1.0) - 

S*PI/LY)**2.0)**2.0)/(ROA*ACHA) 

_  KXM*( (2. 0*EM+1 . 0) *PI/LX)  —  -  - 

KYN-( (2.0'ENI 1 .0)«PI/LY) 

C  COEFICIENT  CALCULATES  THE  INPUT  ARRAY  TO  THE  FFT  (COEF)..  - - - 

CALL  COEFICIENT  (N.HKX.HKY.KXl ,  KYI ,  PI , KXM.KYN, 

_ &  K ,  KC ,  KPF , Z , KX , KY , LX , LY , EM , EN , COEF )  - . . 

C  FFT2B  IS  AN  IMSL  SUBROUTINE  WHICH  CALCULATES  THE  2~D  FFT  OF 

C  A  SET  OF  FOURIER  COEFFICIENTS  . - . -  - 

C  COEF( I ,  J)  IS  THE  INPUT  ARRAY 

C_A(I,J)  IS  THE  OUTPUT  ARRAY  -  --  - 

CALL  F2T2B  ( NR COE F , NCCOEF , A , LDA , COEF , LDCOEF , WFFl , WFF2 , 

. .  HWC.CPY)  ...  . -  —  -  -  : 

DO  400  I-1,N 

_  DO  400  J-l.N  -  - -  - - 

B(  I ,  J)-A( I ,J) 

_  .  _  .  B1(I,J)-(1.0/(2.0*PI)**2.0)*HKX*UKY*( (“1)**I)*( (-1)*?J)  .  - 

S«B(I,J)*(COS((2. 0*EM+1 . 0 ) «PI *XO/LX ) »COS( ( 2 . 0*EN+1 .0) 

_  4*PI*YQ/LY)/(0M£GAMNSQ-QM£GA**2. 0) )  -  - 

BT(  I , J)-BT( I , J)+B1( I , J) 

400  _  CONTINUE  - 

200  CONTINUE 

DO  500  I-l.N  a  .  . . -  - 

DO  500  J-l.M^® 

. .  WB(I.J)r((->r0)«F*RO*QMEGA»*2.0*bT(I,J))/{MP«DB)  -  - 

WBREAL(I ,J)-WB{ I ,J) 

C  WBREAL(I.J)  IS  THE  REAL  PART  OF  THE  PLATE-B  RESPONSE  - 

C  WB(I,J)  IS  THE  COMPLEX  PLATE-B  RESPONSE 

500  CONTINUE 

THE  FOLLOWING  TWO  DO-LOOPS  CONVERT  THE  OUTPUT  PLATE  RESPONSE 
IN  THE  X-Y  PLANE  SO  THAT  THE  RESPONSE  SNOWS  THE  CENTER  OF  THE 
PLATE  IN  THE  CENTER  OF  THE  ARRAY 

1  2  I  I  2 
3  4  |  3  4 
X - 

.12  112  - 

3  4|34 

Y  -  •  - 

DO  600  I-l.N 
DO  600  J-l.N/2 
TEMP-WBREAL  ( I ,  J ) 

WBREAL ( I . J ) • WBREAL ( I . N/2 ♦ J ) 

WDREALf  I ,  N/2U )  -TEMP 
600  CONTINUE  _ 


DO  700  J-1,N 

_ DO_  700  lrlM/2 _ 

TEMP«WBREAL( 1 , J) 

_ HBR£AL(I,  J)?WBREAL(N/2tX,J) _ 

WBR£AL(N/2*I,J)-TEMP 


AMPB-O.Q 

_ CDUNTrO.O _ _ _ _  _ 

C  CALCULATE  THE  AMPLITUDE  OF  VIBRATION  IN  THE  PORTION  OF  THE  PLATE 

C  -  SHADOWED  BY  THE  UPPER  PLATE _ . _ _ _ _ _ 

DO  800  I-(N/2)-(KXl»LX/(4.0«PI)),(N/2)+(KXl*LX/(4.0»PI)) 

- DQ_fl0Q_Jr(N/2)-(KYl;LY/(4 .0?E1) ) , (N/2 )+(KYl?LY/(4*Q?PI) 1 _ 

AMPB-AMPB+WBREAL ( I , J) »*2 . 0 

_ COUNT=OOUNI*l-  a _ _ _ 

800  CONTINUE 

C — AMPB  IS  THE  SUM  OE  THE  SQUARES  OF  THE  .PLATE -RESPONSE  OF  THE  AREA 
C  OF  PLATE  B  SHADOWED  BY  PLATE  A 

C-..AMPTB  IS  AMPB  DIVIDED  BY  THE  NUMBER  OF.  POINTS _ _ 

AMPTB-AMPB/COUNT 

C — CALCULATE— THE  -VIBRATION  RESPONSE  OF.  PLATE -A _ _ 

DO  250  EM-NLO,NUP 

_ DO. 250  ENiNULNUP. _ _ _ _  _ 

OMEGAMNSQ-DSTAR*  ( ( ( (2 . 0«EM+1 . 0 ) «PI/LX) »«2 . 0+( (2 . 0»EN+1. 0) 

- t*£I/LY)*?2.0)**2.  Q)/{ROA!ACHA)  _ _  _ 

DO  450  I-l ,N 

- D0l_45Q.J-L,N _ _  _ _ 

X-(-LX/2.0)+(LX*(I-l))/N 

- Y-(.=LY/2.  Q)+(LY*(  J~l)  )/N _ _ _ 

A(I,J)-(COS((2.0«EM+1.0)*PI*XO/LX))*(COS((2.0*EN+1.0) 

- S?PL?YD/LY) )* (COS(  (2. Q*EM+1 . 0)*PI*X/LX) )? (CQS( (2 .Q*£N+1.Q) _ 

A*PI»Y/LY) )/(OMEGAMNSO-OMEGA«*2 .0) 

- AT( I , J)rAT( I,I)+A(I , J)  _ _ _ _ 

450  CONTINUE 

.250 _ CONTINUE _ _ 

AMPA-0.0 

_ CQUMteQ.O _ _ _ 

DO  550  I— 1  ,N 

_ DQ—550  J=1,N _ _ 

WA(I,J)-(4.0*F/MP)«AT(I,J) 

- WAREAL(I,J)=WA(I,J)  _  _  _ 

C  WAREAL  IS  THE  REAL  PART  OF  THE  PLATE  RESPONSE  OF  PLATE  A 

C — HA  -IS-IHE .  COMPLEX  PLATE  RESPONSE  OF  PLATE.  A _  _ 

AMP A- AMP A+WAREAL ( I , J) «*2 . 0 

_ COUNEsOOUNT+J— 0 _ _ _ _ 

550  CONTINUE 

_ AMPTA-AMPA/COUNT _ _  _ _ 

C  AMPA  IS  THE  SUM  OF  THE  SQUARES  OF  THE  PLATE  RESPONSE  OVER  PLATE  A 

C_  AMPTA.  IS.  AMPA.  DIVIDED  BY  THE  NUMBER  OE.  POINTS _ _ 

WRITE  (26,321)  AMPTA, AMPTB 

C— OUTPUT  AMPTA  .AND.  AMPTB _ _  _ 

321  FORMAT  (25X, 1  AMPTA-  ',E12.6,'  AMPTB-  *,E12.6) 

C-THE  FOLLOWING  . IF. STATEMENT  TESTS  IF  AMPTB  IS  TOO  SMALL _ 

IF ( AMPTB. LE. { 10. 0**( —90. 0) ) )GO  TO  650 

_ XAU-AHPTB/AMPTA _ _  _ _ _ 

C  CALCULATE  TRANSMISSION  LOSS 

- TL(OMEGAI)-10.0»(ALOC10<1.0/TAU)>_ _  _ 

Xl(OMEGAI ) -FLOAT (OMEGA I ) 

C- OUTPUT  TRANSMISSION  LOSS  _ _ _  _ 

WRITE(26 , 311 )OMEGA/(2 . 0*PI ) ,TL(OM£GAI) 

JU - FORMAT  (F12. 2, 10X,F6 .2) _ _ 

GOTO  1000 

65Q _ TL(OMEGAI  )rl01.  Q _  _ _ _ 

WRITE (25 , 311 )OMEGA/(2. 0*PI ) ,TL(OMEGAI ) 

.1000 _ CONTINUE _  _ _ _ _ 

WRITE(26,411) 

411  P0RMAT(//) 
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END 

SUBROUTINE  COEnCIEOT  ( N .HKX.HKY ,KX1, KYI ,PI , KXM, KYN,  ’ 

_ fcK.KC.KPF,Z,KX,KY,LX,LY,EM.£N,COEF)  _ 

C  COEFICIENT  CALCULATES  THE  FOURIER  COEFICIENTS  FOR  USE  IN  TOE 

C.  FET2B  SUBROUTINE  _ 

REAL  KX(N),KY(N),IMN,HKX,HKY,PI,KX1,KY1, 

_  &KXM . KYN . Z , LX , L Y , K  _ _ 

IMECER  EM,EN,N 

-  COMPLEX  TEMP,KC,KPF,COEF(N,N)  _ 

DO  10  I-1,N 

_ DO  10  J-1,N  _ _  _ 

C  TOE  FFT  USES  TOE  INTERVAL  FROM  -KX 1/2.0  TO  KX1/2.0 

C _ AND  -KYl/2.0  TO  KY1/2.Q  _ _  .. 

KX(I)-( (I-l)*HKX)-(KXl/2.0) 

_  ..  KY(J)-(  (J-l)*IIKY)-(KYl/2.0)  _ 

C  THE  FOLLOWING  IF  STATEMENTS  DETERMINE  1MN 

_  IF  ( (KXM?*2 . 0  .£0.  KX(I)?«2.0) -AND.  _  .  _ 

i(KYN**2.0  .EQ.  KY(J)**2.0) ) 

_ iIMN-LX!LY/4 .0  .  _ ... _ 

IF  ( (KXM**2.0  .EO.  KX(I)**2.0) .AND. 

_ fc(KYN**2.0  .NE.  KY(J)**2.0))  _ 

6 IMN-LX*  KYN*COS ( KY ( J ) *LY/2 . 0 ) • ( ~1 . 0 ) *  *  EN/ 

_  i(KYN**2.0“KY(J)»*2. Q)  _ 

IF  {(KXM*«2.0  .NE.  KX(I)**2.0) .AND. 

_ 4<KYN**2.0  .EQ.  KY(J)!*2.Q) )  .  .  _ 

4lMN-LY*KXM*COS(KX(I)*LX/2.0)*(-1.0)**EM/ 

- 4(KXM«*2.0~KX(I)**2.0)  . _ _ 

IF  ({KXM**2.0  .NE.  KX(I)»*2.0) .AND. 

- 4(XYN**2.0  .NE.  KY(J)**2.0) )  ... _ .. 

4IMN-4 . 0*KXM«KYN*COS<  XX{ I ) *LX/2 .  0 )  • 

- 4COS(KY(J)*LY/2.0)*(-L.Q)**£M*(-l.a)**EN/  _  ... 

4((KXM**2.0-KX(I)**2.0)*(KYN»*2.0- 

—  4KY(J)**2.Q))  .  . . _ _ 

C  TOE  FOLLOWING  IF -STATEMENT  TESTS  IF  THE  EXPRESSION  IS 

C_  NEGATIVE  FOR  USE  IN  THE  COMPLEX  EXXPRESSION  FOR  COEF(I.J)  _ 

IF  (K**2.0-KX(I)**2.0-KY(J) 

_ 4**2. 0) .  55.  .65,  65  ._...  .  ...  _ 

55  TEMP-CMPLX(((ABS(K**2.0-KX(I)*«2.0- 

_ 4KY(J)**2.Q))**.5)*(-Z),Q.Q)  _ _  _ 

GOTO  75 

65 _ TEMP-CMPLX(Q.  0, ( (K**2.0~KX(  I) **2. 0~  _ 

4KY(J)**2.0)**0.5)*Z) 

.75  .  CONTINUE  ...  _  _ 

C  CALCULATE  COEF 

_  COEF(I,J)-CEXP(TEMP)*IMN/  _  _ 

4( (KC**2.0-KX(I)**2.0-KY(J)»*2.0)«» . 5* 

_  4((KX(I)**2.0+KY(J)**2.0)»«2.0-KPF))  . 

10  CONTINUE 

_ _ RETURN  ...  .. 

END 

C...  . 

SUBROUTINE  OMEGAMNSUB ( NIO , NUP , OMEGA . ROA , ACH A , DA , LX , LY ) 

C_  OMEGAMNSUB  COMPUTES  TOE  UPPER  AND  LOWER  MODES  WHICH  THE  PLATE 
C  RESPONSE  IS  SUMMED  OVER.  NLO  IS  THE  LOWER  LIMIT  AND  NUP  IS  TOE 
C_  UPPER  LIMIT. 

INTEGER  NLO, NUP, Nil 

- -  REAL  OMEGA , ROA , ACHA , DA , LX , LY ,  OMEGAMNO , OMEGAMN1 , QMEGAMN2 , PI 

PI-2. 0*ASIN( 1.0) 

C _ INITIALIZE  AND  INCREMENT  NUB  _  .  _ 

NLO— 1 

..  15fl  ..  NLO-NLO+1  .... _  _ _ 

C  CALCULATE  RESONANT  OMEGA  FOR  TOE  EN-NIO,  EM -NLO  MODE 

- OMEGAMNO- (DA/ (ROA* ACHA) ) •• . 5*( ( (2. 0*NLO+1. 0) *PI/LX)**2. 0 

**((2.0*NUH1.0)*PI/LY)**2.0) 

C  IF  THIS  RESONANT  OMEGA  IS  LARGER  THAN  OMEGA  GO  ON  TO  258 


C  IF  NOT  INCREMENT  NLO  AND  REPEAT 

- IE_(QMEGAHNQ.GT.QMEGA)  GOTO  258 _  _ 

GOTO  158 

C _ INITIALIZE -ANUJNCREMENIL  NUB _ _ 

258  NUP— 1 

.358 _ NUPsNUP+1 _ _ _ 

C  CALCULATE  THE  RESONANT  OMEGAS  FOR  EN-NUP,  EM-0  AND  EN-0?~EM-NUP 

C _ MODES.. AND JTEST  IF  LARGER  THAN  OMEGA _  _ _  _ _ 

OM£GAMN1-(DA/(ROA*ACHA) )«» . 5»( ( (2. 0»NUP+1.0)«PI/LX)»*2.0 

_ i+(PI/LT)*!2.0) _ _ _ _  _ 

OM£GAMN2-(DA/(ROA*ACHA) )•• . 5* ( { PI/LX)* *2. 0 

- &+(  (2.Q;NUP.+1.0)*PI/LY);?2.0)  _ _ _  _ 

C  IF  LARGER  THAN  OMEGA  GO  ON  TO  458  IF  NOT  INCREMENT  NUP  ,  REPEAT 

- IF(  (QMEGAMN1.GT. OMEGA)  .OR.  (QMEGAMN2.CT. OMEGA) )  GOTO  458 _ 

GOTO  358 

458  _ CONTINUE _ _  _ _ _ _ 

C  ADD  THREE  MODES  TO  UPPER  LIMIT 

_ NUP=NUE*3 _ 

Nll-NLO 

C — SUBTRACT  THREE .  MODES.  FROM.  THE  .LOWER.  LIMIT,  7.FRQ  IS  THE  LQHESTLUL 
NL0-NLO3 

_ IF  (Nil.  EX).  0)  NLQ-Q _ _  _ 

IF(Nll.EQ.l)  NLO-O 

- IE(N11.E0.2)  NLO-0 _ _ _ _  _  _ 

WRITE  (6,558)  NLO,NUP 

.558 _ FORMAT.  (215) _ _ _  _ 

RETURN 


n  o 


PROGRAM  PRTL1000 

CL.  DETERMINES  THE  TRANSMISSION  LOSS  BETWEEN  FLAXE-A  AN11  PLAIE-B-* 
r  USING  RMS  VIBRATION  LEVELS  OF'  EACH  PLATE 

WRITTEN  BY  MICHAEL  F.  SHAW 


INTEGER  OMEGA  I ,  OMEGAB,  OMEGAT,  INC,  EM,  EN,  I,  J, 

_ iNRCOEF , NCCOEF , LOCOEF , LDA , N , NLQ , NUP , Nil  .  _ 

PARAMETER  (N-128) 

_ COMPLEX  .  DSTAR , QMEGAMNSQ , COEF ( N , N )  , A(N,N)  ,B1(N,N) _ _  - 

tBT(N,N) , KPF, KC,WA(N,N) ,WB(N,N),AT(N,N),B(N,N), 

_  AWC(N.N)  _ 

REAL  WFFl(4«N+15),WFF2(4«N+15) ,CPY(N,N) 

REAL  WBREAL(N.N) ,WAREAL(N,N) , WORK! (4000) ,WORK2(100) . 
i AMP A , AMPB , TAU , TL ( 200 ) , COUNT , OMEGAl , OMEGA2 , XI (200) , OMEGART 
REAL  DA, DB. ETAA, ETAB, LX. LY. ROA. R013, ACHA, ACUB,KX(N),KY(N), 

&KX1 , KYl , KXM, KYN, F, RO ,MP , XO , TO , IMN, Z , PRA , GAMMA, MU , NU,C,E, 

_ &tlKX .  HKY ,  PI ,  K ,  KP ,  OMEGA ,  QMEGAMN0 ,  QMEGAMN1 ,  OMEGAMN2  -  -  - 

INTRINSIC  CMP LX 

_  EXTERNAL  F2T2B.FFTCX  _  _ 

C  TOE  NEXT  TWO  LINES  ALLOW  N  TO  BE  LARGER  THAN  64  POINTS 

_  COMMON  /WORKSP/  RWKSP  .  -  - 

REAL  RWKSP (8884) 

C  . . VARIABLE  LIST 

C  DA~-  BENDING  llIGIDITY  OF  THE  FINITE  PLATE 

C _ DB  -BENDING  RIGIDITY  OF  THE  INFINITE  PLATE _ : - 

C  ETAA  -  DAMPING  COEFICIENT  FOR  TOE  FINITE  PLATE 

C.  _  ETAB .r  DAMPING  COEFICIENT  FOR  THE  INFINITE  PLATE  _ 

C  LX  t  LY  -  DIMENSIONS  OF  TOE  FINITE  PLATE 

C _ ROA  -  DENSITY  OF  THE  FINITE  PLATE  _ 

C  ROB  -  DENSITY  OF  TOE  INFINITE  PLATE 

C  ACHA  r  THICKNESS  OF  THE  FINITE  PLATE  _ 

C  ACHB  -  THICKNESS  OF  THE  INFINITE  PLATE 

C _ RQ  =  DENSITY  OF  FLUID  _ 

C  XO  *  YO  -  TOE  POINT  WERE  TOE  INPUT  FORCE  F  IS  APPLIED 

C _ Z  -  DISTANCE  BETWEEN  THE  PLATES  _ _ 

C  PRA  -  PRANDEL  NUMBER  FOR  FLUID 

C.  _  GAMMA  -  SPECIFIC  HEAT  RATIO  FOR  FLUID  _ _ 

C  MU  -  VISCOSITY  OF  FLUID 

C  _  NU  -  POISON’S  RATIO  .  . .  . . 

C  C  -  SPEED  OF  SOUND  IN  FLUID 

C  EA  r  SHEER  MODULUS  OF  THE  FINITE  PLATE  _ 

C  EB  -  SHEER  MODULUS  OF  THE  INFINITE  PLATE 

C  _KX1  -  UPPER  WAVENUMBER  LIMIT  IN  FINITE  PLATE  IN  X  DIR _ 

C  KYl  -  UPPER  WAVENUMBER  LIMIT  IN  FINITE  P1ATE  IN  Y  DIR 

PARAMETER  (F- 1.6) 

C  .  INPUT  FORCE  TO  THE  FINITE  PLATE  IS  ASSUMED  TO  BE  A  UNIT  FORCE  (1)  — 

OPEN  (UNIT-25,  FILE-’ PR1000.DAT’ ,  STATUS- ’OLD’ ) 

OPEN  (UNIT-26,  FILE? ’PRTL1000. OUT’ ,  STATUS- ’ NEW* )  - -  - 

C  ANOTHER  LINE  SO  N  CAN  BE  LARGER  THAN  64 
_  .  CALL  IWKIN( 8084) 

C  READ  INPUT  DATA  FROM  PR. DAT 

READ  (25,15)  LX. LY.XO.YO, ACHA. ACHB  .  . . - 

15  PORMAT(4F8. 3.2F8.6) 

READ  (25,25)  ROA, BOB, RO. ETAA. ETAB, GAMMA. PRA  _  - 

25  FORMAT (2F8. 2, F8. 3.2F8. 7 ,2FH. 6 ) 

_  READ  (25.35)  MU.C.EA.EB.KXl.KYl  -  - 

35  PORMAT(F8.7.F8.3,2E8.3.2F8.2) 

_  READ  (25,45)  NU. OMEGAl, OMEGA2. I NC.X, OMEGA  -  ..  - 

4b  FORMAT (F8. 4 ,F8 . 1 ,F16 . 1 , I8,F8. 3 ,F12. 3) 

C  COMPUTE  TOE  VALUE  OF  PI  -  -  - 

PI-2. 0*ASIN( 1.0) 

C  COMPUTE  TOE  BENDING  RIGIDITIES  _ _ 


DA-(EA«ACHA«*3.0)/(12.0«(1.0-NU««2.0)) 

_ DB*?(EB*A01B*;3.  D)/(  12.0* (1. Q-NU£l2.Q)  ) _ 

DSTAR-CMPLX < DA, DA* ETAA) 

_ M£-ROA*ACHA*LX*LY  _ _  _ _ _ 

WRITE(26,11)  LX,LY 

JJ _ FORMAT}  L  PARAMETERS  OF  FINITE  PLATE',/,  i.  DIMENSIONS  0FPLATE1 

i,‘  (IX)  1 , F8 . 3 , ‘  (LY)  ' , F8 . 3 , '  METEnS') 

_ _ HRITE<26,21)  XQ,YO  —  _..  __ _ 

21  FORMAT} 1  DRIVING  POINT  XO  ' ,F8.3,'  YO  \F8.3,' 

_ &  METERS') _  _  —  ... _  .  - 

WRITE<26,31)  F 

JJ _ FORMAT}  1  DRIVING  FORCE .  FB .  3 , 1-  N ' ) _ _ 

WRITE (26 ,4 1)ACHA 

41 _ FORMAT} '  PLATE  THICKNESS  :  .F8.fi,-*  K* ) _  _ 

WRITE(26,51)  ROA 

J1  _  EORMATC  PLATE  DENSITY  '.F8.3,'.  KG/CUrMETERS  * )  _ _ 

WRITE(26,61 )  EA 

61 _  FORMAT} :  MODULUS .j  OF  ELASTICITY  J ,E8. 3.1  M/SQ-METERS * ) _ 

WRITE} 26, 71)  DA 

2J _ FORMAT} BENDING  RIGIDITY.  ...!., E12 .4 )  —  _ _ _ 

WRITE} 26 ,81)  ETAA 

flj _ FORMAT f. .DAMPING  COEFFICIENT F8 . 4 )  _ . .. _ 

WRITE} 26, 91) 

_91 _ FORMAT}/,  1  PARAMETERS . OF  INFINITE  PLATE',/) . . 

WRITE(26 ,101 )  ACHB 

.101 _ FORMAT} '  PLATE  THICKNESS  ‘  .FB.fi M‘.) - - 

WRITE(26,111)  ROB 

11D__F0RMAT} I_ PLATE  DENSITY-  :,F8. 3,  J _ KG/CURMETERS ' ) _ 

WRITE (26, 121)  EB 

__121 _ FORMAT}  1  MODULUS  OE. ELASTICITY  _'.,EB. 3,  L- N/SQ-METERS L) _ 

WRITE(26,131)  DB 

.131 _ FORMAT}  L  BENDING  RIGIDITY-  '.,E12-4) - . - 

WRITE} 26, 141)  ETAB 

J4D_FORMAT( ' -DAMPING  CQEEICIENT’  ,F8.4) - - 

WRITE (26, 151)  RO 

_151__FQRMAT(/.L  PARAMETEBS  FOR  FLUID'.,/,—  DENSITY — :  ,FB.3, _ 

6  *  KG/CU-METER') 

_ WRITE(2fi,lfil)  MU. - 

161  FORMAT} '  VISCOSITY  \F8.7,'  KG/M-S’) 

_ WRITE} 26. 181)  C -  -  - 

181  FORMAT} '  SPEED  OF  SOUND  ',F8.3,'  M/S') 

_ HRITE(26,191)  PRA - -  - - 

191  FORMAT} '  PRANOTL  NUMBER  ’,F8.3) 

_ WRITE(2fi,201)  GAMMA -  - - 

201  FORMAT}'  SPECIFIC  HEAT  RATIO  \F8.3) 

_ (SITE} 26, 211)  2  . . . .  .  —  -  -  - 

211  FORMAT} '  DISTANCE  BETWEEN  PLATES  ',F8.3,'  M') 

C  .FREQUENCY  IS.  TAKEN  FROM  125  .TO  1000  M2  . - - - 

DO  1000  OMEGAI-21,30 

_ OMEGA =2. 0?PI* (IQ. Q;? (QMEGAI/10. Q) )  -  - 

C  TOE  LIMITS  OF  TOE  FFT.  TOE  FFT  IS  TAKEN  FROM  -KX1/2.0  TO 

C  KX1/2.Q  AND  FROM  -KYI/2.0  TO  KYl/2.0  - 

KX1-201.1 

_ KYI-201.1 _  .  -  ....  - - 

C  TOE  SAMPLE  SPACING  IS  HKX  AND  HKY 
_ _  HKX-JOl/N  .  .  - -  - 

hky-kyi/n 

C  —PARAMETERS  FOB  XFI28  INVERSE  FOURIER  TRANSFORM  SUBROUTINE _ 

NRCOEF-N 

_ NCCOEF-N - - 

LDCOEF-N 

_ 1DA*N - -  - - - 

C  INITIALIZATION  ROUTINES  FOR  TOE  FIT 

_  .  riJJ,  FFTCI(N,WFF1) -  - - - 

CALL  FFTCI (N.WFF2) 

C  K  IS  WE  WAVENUMBER  IN  TOE  FLUID _ 


nononnnnno 


K-OMEGA/C 

C-  KC  IS  TOE  COMPLEX  WAVENUMBER  .  _ . _ 

KC-(K)*CMPLX(1.0,  (OMEGA*MU*< (4. 0/3.0)+ ( (GAMMA-1. 0)/PRA) 

_  i)/(2.0*RO»C?!2.0)))  .  .  _ 

C  KP  IS  TOE  PLATE  WAVENUMBER 

_  KP“ ( ROB* ACHB*QM£GA*  *2 . Q/DB ) * ’ . 25  _ 

C  KPF  IS  TOE  COMPLEX  PLATE  WAVENUMBER* *4 . 0 

KPF-(ROB*ACHB*CtM£GA**2.0/DB)?CMPLX(i.0,ETAB) _ 

00  100  I-1,N 

_  DO  100  Jrl,N  _  _ 

BT(I,J)-0.0 

_  AT(I,J)-0.0  .  .  _  -  - 

100  CONTINUE 

C  QMEGAMNSUB  CALCULATES  TOE  UPPER  AND  LOWER  MODES  NUP,  NLO  -  - 

CALL  OMEGAMNSUB  ( NLO ,  NUP ,  OMEGA ,  ROA .  ACI I A ,  DA ,  LX ,  LX ) 

C  TOE  PLATE  MODES  ARE  TAKEN  FROM  EM  AND  EN  -  NLO  TO  NUP  _.  . 

DO  200  EM-NLO.NUP 

DO  200  EN-NLO.NUP  - - 

C  OMEGANMSQ  IS  TOE  RESONANT  OMEGA**2 

_  OMEGAMNSQ-DSTAR* ( ( ( (2 . 0*EMH.  0)  *PI/LX)  ??2.  Q+(  (2.0?ENtl.O) - 

S*PI/LY)**2.0)**2.0)/(ROA*ACMA) 

_  KXM-((2.0*EM+1.0)*PI/LX)  - 

KYN“( (2. 0*EN+1 . 0) *PI/LY) 

C..  COEFICIENT  CALCULATES  TOE  INPUT  ARRAY  TO  TOE  FFT.(COEF) - 

CALL  OOEFICIENT  (N,HKX,HKY,KX1 ,KYl ,PI , KXM,KYN, 

_  4K.KC.KPF,Z.KX,KY.LX.LY,EM.EN,C0EF)  - 

C  FFT2B  IS  AN  IMSL  SUBROUTINE  WHICH  CALCULATES  THE  2-D  FFT  OF 

C_  A  SET  OF  FOURIER  COEFFICIENTS  .  - : - 

C  COEF( I , J)  IS  TOE  INPUT  ARRAY 

C  A( I , J )  IS  TOE  OUTPUT  ARRAY  .  _ 

CALL  F2T2B  ( NRCOEF , NCCOEF , A , LDA , COEF , LDCOEF , WFFl , WFF2 , 

_ AWC/CPY) .  ....  ... - 

DO  400  I-1,N 

DO  400  J-1,N  ..  - - 

B(I,J)-A(I,J) 

_ B1(I.J)-(1.0/(2.0«PI)**2.0)*HKX*HKY*((-1)«*1)*((=1);*J.) _ 

4*B(I,J)*(COS((2.0*EM+1.0)*PI*XO/LX)*COS((2.0*EN+1.0) 

_ t*Pl?YO/LY)/(QM£GAMNSO-OM£GA**2. 0) )  - - 

BT(I,J)-BT(I,J)+B1(I,J) 

.400  CONTINUE  _  - 

200  OONTINUE 

I.  DO  500  Isl.N  _  _  ...  .  - 

DO  500 

.  _  .  WB(I,J)-((-\fD)?F*RO?QMEGA**2.0*BT(I,J))/(MP»DB) - 

WBREAL(I,J)-WB(I,J) 

C_  WBREAL(I.J)  IS  TOE  REAL  PART  OF  THE  PLATE-B  RESPONSE _ 

C  WB( I , J)  IS  TOE  COMPLEX  PLATE-B  RESPONSE 

!  500  CONTINUE  -  -  - 

TOE  FOLLOWING  TWO  DO-LOOPS  CONVERT  TOE  OUTPUT  PLATE  RESPONSE 

IN  TOE  X-Y  PLANE  SO  TOAT  TOE  RESPONSE  SHOWS  TOE  CENTER  OF  TOE. - 

PLATE  IN  TOE  CENTER  OF  TOE  ARRAY 

12  112 
3  4  ]  3  4 

X  - 

1  2  I  1  2  ..... 

3  4  I  3  4 

_  Y  •  -  - - 

DO  600  I-1,N 

DO  600  J-l.N/2  - - —  - - 

TEMP  *408  CAL  ( I  ,J) 

WBREAL(  I .  J )  *WDR£AL(  I  .N/2+J }  -  ...  . 

WBHEAL(I ,N/2+J)-TEMP 

600  CONTINUE  - - 

DO  700  J-l.N 
DO  700  I-l.N/2 


TEMP-WBREAL(I,J) 

_  WBREAL(I,J)-WBREAL(N/2+I, J) _ 

WBREAL(N/2+I,J)-TEMP 

700 _ CONTINUE _ —  .  _ _ _ 

AMPB-0.0 

_ _ COUNT-O.Q . .  ...  —  - 

C  CALCULATE  THE  AMPLITUDE  OF  VIBRATION  IN  TOE  PORTION  OF  TOE  PLATE 

C _ SHADOWED  BY  TOE  UPPER  PLATE  ....  .  _  .  - - - 

DO  800  I-(N/2)-(KXl»LX/(4.0«PI)),(N/2)+(lCXl*LX/(4.0*PI)) 

_ DO  800 . J- ( N/2 ) - ( KY1*LY/ <  4 . 0«  PI ) ) , ( N/2 )  +  ( KYI *LY/ ( 4 . Q»PI  _ 

AMPB- AMPB4WBREAL ( I ,  J  >  «  •  2 . 0 

_ COUNT-CQUNT+l.O _  _..  - - 

800  CONTINUE 

C  AMPB  IS  THE  SUM  OF  TOE  SQUARES  OF  TOE  PLATE  RESPONSE  OF  TOE  .AREA 
C  OF  PLATE  B  SHADOWED  BY  PLATE  A 

C  AMETB  IS  AMPB  DIVIDED  BY  TOE  NUMBER  OF  POINTS _  _ _ 

AMPTB-AMPB/COUNT 

C_  CALCULATE  TOE  Y1 BRAT ION  RESPONSE  OF  PLATE  A -  - 

DO  250  EM-NLO , NUP 

_ DQ-25Q  ENsNLO.NUP  -  _  -  - 

OMEGAMNSQ-DSTAR« ( ( ( ( 2 . 0*EM+1 . 0 ) *PI/LX) *«2 . 0+ ( ( 2 . 0*EN+1 . 0) 

_ i«EI/LY)**2.Q)  5!2. 0)/(ROA*ACHA)  - -  ...  - 

DO  450  I-l.N 

_ DQ—45Q- J“1.N _  _ _  _ 

X-(-LX/2.0)+(LX»{I-l))/N 

_ Y=(rl.Y/2.Q)t(LY?(J-l))/N . . .  . 

A(I,J)-(COS((2.0»EM+1.0)*PI*XO/LX))*(OOS({2.0*EN+1.0) 

_ i*PHYO/LY).)*<OOS{  (2.0*EM+1.0)*PI*X/LX)  )*(COS(  (2.0*EN+1.0) - 

4*PI«Y/LY) )/(OMEGAMNSQ-OMEGA«*2 .0) 

_ AT(I*J)sAT(L,J)+A(I,J)  ... - - 

450  CONTINUE 

250  CONTINUE - -  - -  - 

AMPA-0.0 

_ DQUNfc.0 . 0 . . . -  - - - 

DO  550  I-l.N 

_ EO..S50  J-1,N _  —  -  - - - 

WA( I , J)-(4 . 0*F/MP) *AT( I ,  J) 

_ HAREAL(I.J)=HA(I,J)  -  - -  - 

C  W AREAL  IS  TOE  REAL  PART  OF  TOE  PLATE  RESPONSE  OF  PLATE  A 

C _ WA  IS  TOE-COMPLEX  PLATE  RESPONSE  OF  PLATE  A  -  - 

AMPA-AMPA+WAREAL( I , J) ••2 . 0 

_ CQUNTrCOUNTtl.D _  .  —  - -  - - 

550  CONTINUE 

_ AMPTA-AMPA/COUNT _  _  _ _  _ 

C  AMPA  IS  TOE  SUM  OF  TOE  SQUARES  OF  TOE  PLATE  RESPONSE  OVER  PLATE  A 

ir  ampta  is  AMPA  DIVIDED  BY  TOE  NUMBER  OF  POINTS  ..  _ 

WRITE  (26,321)  AMPTA, AMPTB 

C _ OUTPUT  AMPTA- AND  AMPTB  .  _  ..  . - - - - 

321  FORMAT  (25X,*  AMPTA-  \E12.6,*  AMPTB-  \E12.6) 

C  THE  FOLLOWING  IF  STATEMENT  TESTS  IF  AMPTB  IS  TOO  SMALL  _ 

IF( AMPTB. LE. (10.0**(-90.0) ) )GO  TO  650 

_ _ TAU-AMPTB/AMPTA  ..  - 

C  CALCULATE  TRANSMISSION  LOSS 

_  TL ( OMEGAX ) “10 . 0* ( ALOG1 0 { 1 . 0/TAU ) )  .  _  - 

XI (OMEGA I )-FL3AT(0MEGAI ) 

C  OUTPUT  TRANSMISSION  LOSS  - - -  - 

WRITE(26 , 311 )OMEGA/( 2 . 0«PI ) ,TL(OMEGAI ) 

J1J _  FORMAT(Fl2.2, 10X,F6 .2)  - -  - 

GOTO  1000 

650 _ TL(OM£CAI)-101.0  .  - - 

WRITE ( 26 , 311 )OHEGA/ ( 2 . 0«PI ) ,TL(OMEGAI ) 

JL000 _ CONTINUE.  .  .  _  .  -  - - 

END 


SUWOUTXNE  COEFieXENT  (N,HKX,HKY,KX1,KY1,PI,KXM,KYN, 
*K, KC , XPF , 2 , KX , KY , LX , LY , EM , EN . COEF) _ 


C  COEFICIENT  CALCULATES  TOE  FOURIER  COEFICIENTS  FOR  USE  IN  TOE 

C _ FFT2B  SUBROUTINE  .  .  ...  _ _ _ 

REAL  KX(N),KY(N),IMN,HKX,HKY,PI,KX1,KY1, 

_ 4KXM,KYN.Z,LX,LY,K  .  _ 

INTEGER  £M , EN , N 

_ COMPLEX  TEMP, KC,KPF,COEF(N,N)  - _ 

DO  10  I-1,N 

_  DO  10  J-l.N  _  _ 

C  THE  FFT  USES  THE  INTERVAL  FROM  -KX1/2.0  TO  KXl/2.0 

C  AND  -KYI/2.0  TO  KYI/2.0  .  _ 

KX<I)-<(I-l)*HKX)-<KXl/2.0) 

_ KY( J)-( (J-l) *HKY)“(KYl/2. 0)  - 

C  TOE  FOLLOWING  IF  STATEMENTS  DETERMINE  IMN 

IF  ( (KXM**2.0  .EO.  KX(I)?«2.Q) . AND.  _  - - 

4(KYN**2.0  .EQ.  JCY(J)«*2.0)) 

4lMN-LX*LY/4 . 0  ....  - -  - - 

IF  ( (KXM**2.0  .EO.  KX(I)«*2.0) .AND. 

_  MKYNV2.0  ,NE.  KY(J)»*2.0>)  - 

i IMN-LX  *KYN*COS ( KY ( J ) • LY/2 . 0 ) • ( -1 . 0 ) « • EN/ 

i(KYN**2.0-KY{ J)**2.0)  ..  - - 

IF  ( (KXM**2.0  .NE.  KX(I)*»2.0) .AND. 

_ fc(KYN**2.0  .EQ.  KY(J)“*2.Q) )  ..  .  - 

4IMN-LY* KXM*COS ( KX ( I ) * LX/2 . 0 ) • { -1 . 0 ) • • EM/ 

_ i(KXM*?2.0-KX(I)**2.0)  .  .  - 

IF  ( (KXM**2.0  .NE.  KX(I)**2.0) .AND. 

_ &(KYN**2.0  .NE.  KY(J)**2.0))  - 

iIMN-4 . 0*KXM«KYN*COS ( KX ( I ) »LX/2 . 0) • 

_  iCOS(KY(J)?LY/2-0)*(-1.0)**EM*(rl.0)**EN/  - 

& ( ( KXM*  *  2 . 0-KX ( I ) *  *  2 . 0 ) • (  KYN" *2 . 0- 

_ m(J)?*2.0))  ...  ...  - 

C  TOE  FOLLOWING  I /-STATEMENT  TESTS  IF  TOE  EXPRESSION  IS 

C_  NEGATIVE  FOR  USE  IN  TOE  COMPLEX  EXXPRESSION  FOR  COEEfl^J) - 

IF  (K«*2.0-KX(I)**2.0-KY(J) 

_ t**2.0)  55,  65,  65  -  - 

55  TEMP-CMPLX(((ABS(K**2.0-KX(I)«*2.0- 

_ 6XY( J)**2.0) )•• .  5)*(~Z) ,0.0)  _  - 

GOTO  75 

_£5 _ T£MPj<MPLX<Q.0.((X?*2.0-KX(I)**2.0-  - 

iKY(J)**2.0)**0.5)*Z) 

_75 _  CONTINUE  - 

C  CALCULATE  COEF 

_  COEF(I.J)-CEXP( TEMP)* IMN/  - 

i( (KC**2 . 0-KX(I) **2. 0-KY( J)**2 . 0) ** . 5* 

_  &( (KX(I)?*2.0+KY(J)**2.0)**2.0-KPF) )  - 

10  CONTINUE 

_ RETURN  -  - - 

END 

C _  _ 

SUBROUTINE  OMEGAMNSUB  ( NLO ,  NUP ,  OMEGA ,  ROA ,  ACHA ,  DA ,  LX ,  LY ) 

C_  QMEGAMNSUB  COMPUTES  THE  UPPER  AND  LOWER  MODES  VWICH  TOE  PLATE - 

C  RESPONSE  IS  SUMMED  OVER.  NLO  IS  TOE  LOWER  LIMIT  AND  NUP  IS  TOE 
|C  UPPER  LIMIT. 

INTEGER  NLO, NUP. Nil 

REAL  OMEGA.ROA, AQlA.DA,LX.LY,QMEGAMN0 ,QMEGAMNl,QMEGAMN2 ,PX  —  - 

PI-2. 0*ASIN< 1.0) 

C  INITIALIZE  AND  INCREMENT  NUP  -  - 

NLO— 1 

158  NLO-NLO+1  - - 

C  CALCULATE  RESONANT  OMEGA  FOR  TOE  EN-NIO.  EM -NLO  MODE 

QMECAMN0-(DA/(ROA*ACHA) ) •* .5* ( < (2 . O'NLOtl. 0) *PI/LX) ?*2. 0 -  _.  - 

*♦( (2. 0*NLO*1 .0)*PI/LY)**2 . 0) 

C  IF  THIS  RESONANT  OMEGA  IS  LARGER  TOAN  OMEGA  GO  ON  ID  258 - - 

C  IF  NOT  INCREMENT  NLO  AND  REPEAT 

IF  (OMEGAMNO.CT. OMEGA)  GOTO  258  _  - - 

GOTO  158 

C  INITIALIZE  AND  INCREMENT  NUP  _ 


258  NUP— 1 

.258 _ NUE-NUEJJ _ _ _ 

C  CALCULATE  THE  RESONANT  OMEGAS  FOR  EN-NUP,  EM-0  AND  EN-O,  EM-NUP 

C— MODES -AND  JEST  IF  LARGES  TUAN  OMEGA ....  .  _  _ 

OMEG AMN 1 - ( DA/ ( RO A « ACH A ) ) • • . 5 • < ( { 2 . 0 • NUP ♦ 1 . 0 ) • P I /LA ) • • 2 . 0 

_ it(PI/LY)»;2.0) _ ...  ..  . . .  . 

OM£GAMN2-(DA/(ROA«ACHA) ) •« . 5*( (PI/LX)««2. 0 

_ i+((2.0!NUE+l.Q)!PI/LY)*«2.0)  _  _  _ 

C  IF  LARGER  THAN  OMEGA  GO  ON  TO  458  IF  NOT  INCREMENT  NUP  ,  REPEAT 

_ IF(  (QMEGAMN1.GT. OMEGA)  .OR.  (QMEGAMN2.GT. OMEGA) )  GOTO  458 _ 

GOTO  358 

_158 _ CONTINUE _ _ _ 

C  ADD  THREE  MODES  TO  UPPER  LIMIT 

_ NUE-NUP+3 _ _  .  _ _ _ _  _ 

Nil-NLO 

C_  SUBTRACT  THREE  MODES  FROM  THE  LOWER  LIMIT,  ZERO  IS  UtE  LOWEST.N) 
NLO-NLO-3 

_ IF(Nll.EO.O)  NLOrO _ _ _  _ _ 

IF(Nll.EQ.l)  NLO-O 

_ IF'(N1L.EQ.2)  NLO-O _ _ _ _ _  _ 

WRITE  (6,558)  NIC, NUP 

_558 _ EORMAT  (215) _ _  _ 

RETURN 


on  i  oi  i  i  :  ,  noon 
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PROGRAM  PRTL100 

C _ DETERMINES  TOE  TRANSMISSION  LOSS  BETWEEN  PLATEtA  AJJn  PLATE-B  . _ 

USING  RMS  VIBRATION  LEVELS  OF  EACH  PLATE 

""WRITTEN  BY  MICHAEL  F.  SHAW 

’  INTEGER  OMEGAI,  OMEGAB,  OMEGAT,  INC,  EM,EN^I ,  j, 

_ 4NRCOEF , NCCOEF , LDCOEF , LDA , N ,  NLO , NUP , Ni 1  _ . _ 

PARAMETER  (N-64) 

_ COMPLEX  DSTAR,QMEGAMNSQ,COEF(N,N) ,A(N,N) ,B1(N,N), _ 

4BT(N,N) ,KPF,KC,WA(N,N) ,WB(N,N),AT(N,N),B(N,N), 

_ 4WC(N,N)  ...  _ 

REAL  WFF1(4*N+15) ,WFF2(4«N+15) ,CPY(N,N) 

_ REAL  WBREAL(N.N) , WAREAL(N.N) , W0RK1(  4000)  ,WORK2(100) ,  -  - 

&  AMP  A ,  AMPB ,  TAU ,  TL  ( 2  00 ) ,  COUNT ,  OMECAl ,  OMEGA2 ,  XI  (200) ,  OMEGART 

_  REAL  DA.DB,ETAA,ETAB,LX.LY,ROA.ROB,ACHA.ACHB.KX(N),KY(N), 

4KXl,KYl,KXM,KYN,F,R0,MP,X0, YO,IMN,Z,PRA. GAMMA, MU, NU,C,E, 

..  ..  4HKX ,  HKY .  PI ,  K ,  KP ,  OMEGA ,  QMEGAMN0 ,  QMEGAMN1 ,  OMEGAHN2 
INTRINSIC  CMP LX 

_  EXTERNAL  F2T2B , FFTCI  -  _ _ _ - 

THE  NEXT  TWO  LINES  ALLOW  N  TO  BE  LARGER  THAN  64  POINTS 

COMMON  /VDRKSP/  RWKSP  ...  .....  _ _ _ 

REAL  RWKSP (8884) 

. VARIABLE  LIST  ~  .  . .  "" 

C _  _  _ 

C  DA  -  BENDING  RIGIDITY  OF  TOE  FINITE  PLATE 

C _ DB  -BENDING  RIGIDITY  OF  TOE  INFINITE  PLATE  _ ^ _ 

C  ETAA  -  DAMPING  OOEFICIE^T  FOR  TOE  FINITE  PLATE 

C ETAB  r  DAMPING  COEFICII.T  FOR  TOE  INFINITE  PLATE _ _ 

C  LX  4  LY  -  DIMENSIONS  Of  TOE  FINITE  PLATE 

C _ ROA  r  DENSITY  OF  TOE  FINITE  PLATE  _  . 

C  ROB  -  DENSITY  OF  TOE  INFINITE  PLATE 

C _ ACHA  -  THICKNESS  OF  THE  FINITE  PLATE  .  _ _ _ 

C  ACHB  -  THICKNESS  OF  TOE  INFINITE  PLATE 

C_  RQ  -  DENSITY  OF  FLUID  _ _ _ _ 

C  XO  A  YO  -  TOE  POINT  WERE  TOE  INPUT  FORCE  F  IS  APPLIED 

C  Z  DISTANCE  BETWEEN  TOE  PLATES  -  - - - - 

C  PRA  -  PRANDEL  NUMBER  FOR  FLUID 

C  _  GAMMA  -  SPECIFIC  HEAT  RATIO  FOR  FLUID  _ _ 

C  MU  -  VISCOSITY  OF  FLUID 

C _  NU_=  POISON'S  RATIO  _  ... _ 

C  C  -  SPEED  OF  SOUND  IN  FLUID 

C EA  -  SHEER  MODULUS  OF  TOE  FINITE  PLATE  _ 

C  EB  -  SHEER  MODULUS  OF  TOE  INFINITE  PLATE 

C_  KX1  -  UPPER  WAVENUMBER  LIMIT  IN  FINITE  PLATE  IN  X  DIR _ 

C  KYI  -  UPPER  WAVENUMBER  LIMIT  IN  FINITE  PLATE  IN  Y  DIR 

PARAMETER  (fal  0) 

C...  INPUT  FORCE  TO  TOE  FINITE  PLATO  IS  ASSUMED  TO  BE  A  UNIT  FORCE  (1) 
OPEN  (UNIT-25,  FILE-' PRI00.DAT' ,  STATUS-' OLD' ) 

_  OPEN  (UNIT-26,  FILE- 'PRTL100.0UT' ,  STATUS- 'NEW'  ) _ 

C  ANOTHER  LINE  SO  N  CAN  BE  LARGER  THAN  64 
CALL  IWXIN<8884) 

C  READ  INPUT  DATA  FROM  PR. DAT 

READ  (25.15)  LX.LY.XO, YO.AdlA.ACUB  _  .  ...  . 

15  FORMAT(4F8. 3,2F8. 6) 

READ  (25,25)  ROA. ROB. RO.CTAA. ETAB, GAMMA, PRA  _ 

25  PORMAT( 2F8 . 2 , F8 . 3 , 2F8 . 7 , 2F8 . 6 ) 

_  ..  READ  (25,35)  MU.C.EA.EB.KXl ,KYl  _ 

35  PORMAT(F6. 7,F8. 3.2E8. 3 ,2F8.2) 

READ  (25,45)  NU.OMEGAl .OMEGA2. INC, Z. OMEGA  -  - 

45  PORMAT(F8. 4 ,F8.1,F16.1,I8,F8.3,F12.3) 

C  COMPUTE  TOE  VALUE  OF  PI  -  .  _  _ 

PI-2. 0*ASIN( 1.0) 

C  COMPUTE  TOE  BENDING  RIGIDITIES 


op 


DA-( EA*ACHA*« 3 . 0 ) /{ 12 . 0« (1 . 0-HU««2 . 0 > ) 

_ -DB?(EB*ACHB*i3.0)/(12.0*(1.0-NU?i2.0)  ) _ _ 

DSTAR-CMPLX  { DA ,  DA*  ETAA) 

_ mp-rqa*acha?lx;ly _  _ 

WRITE(26,11)  LX,LY 

.11 _ PORMATC  PARAMETERS  OF  FINITE  PLATE L  DIMENSIONS  OE-PLAEBl 

k, 1  (LX)  1 ,F8.3, '  (LY)  ‘.Fe.I,'  METERS') 

_ WRITE(26,21)  XQ.YQ .  .  . 

21  FORMAT  ( '  DRIVING  POINT  XO  ' ,F8.3,'  YO  ',F8.3,' 

_ A-METESS')-.. _  _  _  .  .  .  _ 

WRITE(26 , 31)  F 

.  11 _ FORMAT ( •.  DRIVING  FORCE  '_,FB .  3 , '  .M‘  > _  ,  .  _ _ 

WRITE(26,41)ACHA 

41 _ FORMAT ( !  PLATE  THICKNESS  ,FB.6,_i_  M' ) _ _ 

WRITE(26,51)  ROA 

51 _ FORMAT ('  PLATE  DENSITY  * ,  F8 . 3 , KG/CU-METERS ' )  .  _ 

WRITE(26 ,61)  EA 

61 _ FORMAT  (!  MODULUS  OF  ELASTICITY  1  ,E8. 3. N/SQ-METERS1.) _ 

WRITE(26,71)  DA 

J1 _ FORMAT (  -  BENDING  RIGIDITY  ' ,E12.4)  _  _ 

WRITE (26 ,81)  ETAA 

. .81 _ FORMAT ( 1  DAMPING  COEFFICIENT 1  ,FB.4) _  _ 

WRITE(26,91) 

_31 _ FORMAT(/,  1  PARAMETERS  OF  INEINITE  PLATE !_/)_  _ 

WRITE(26 , 101 )  ACHB 

1.01 _ FORMAT ( 1.  PLATE. THICKNESS-  „i, FB.fi,!..  M!  ) _ _ 

WRITE (26, 111)  ROB 

111 _ FORMAT ( 1 -PLATE  DENSITY  ' ,FB.3,i_.  KG/CU-METERS' ) _ 

WRITE (26, 121)  EB 

12) _ FORMAT  ( ’-MODULUS  OF  .ELASTICITY  1.E8.3C N/SO~METERS  *) _ 

WRITE(26,131)  DB 

131  FORMAT ( L.BENDING  -RIGIDITY  '  ,E12. 4) _ _ _ 

WRITE (26 ,141)  ETAB 

-141 _ EQRMAXC . DAMPING  (DEFICIENT  ,F8 .4) . .  . 

WRITE(26,151)  RO 

151 _ FORMAT ( A+l  PARAMETERS  FOR  FLUID.'-/,!  DENSITY.  ',F8.3, _ 

f  KG/CU-METER') 

_ HRITE(26,161)  MU.  . .  .  . . . . 

161  FORMAT ( '  VISCOSITY  \F8.7,'  KG/M-S') 

_ WRITE(26,181)  C  .  .  _ _  _ 

181  PORMATC  SPEED  OF  SOUND  \F8.3,'  M/S') 

_ HR1TE(26 ,191)  PRA.  _ _ _ _ 

191  PORMATC  PRANOTL  NUMBER  \F8.3) 

_ HRITE(26,2QL)  GAMMA  .  _ _  _ 

201  PORMAT( '  SPECIFIC  HEAT  RATIO  ',F8.3) 

_ WRITE (26, 211)  Z_  _  _  _ 

I  211  PORMATC  DISTANCE  BETWEEN  PLATES  ',F8.3,‘  M') 


TOE  FREQUENCY  RANGE  IS  TAKEN  FROM  10  TO  100  Hz 

_ WRITE(26,301) _  _ _ _ _ _ _ 

301  PORMATC  FREQUENCY  (HZ)  TRANSMISSION  LOSS  ') 

_ DQ-1000  DMEGAI-10.20 _ _  _  _ _ _ 

OMEGA-2 . 0»PI • ( 10 . 0»* (OMEGAI/10 . 0 ) ) 

C _ KX1  AND  KYI  —ARE _ _  _ _ 

C  TOE  LIMITS  OF  WE  FFT.  TOE  FFT  IS  TAKEN  FROM  -KX1/2.0  TO 

C_  KX1/2.Q.AND  FROM  -KYI/2.0  TO  KYI/2.0  _ _ 

KX1-62.9 

_ KYl-62.3 . . . . _ 

C  TOE  SAMPLE  SPACING  IS  HKX  AND  HKY 

_ HKX-KXl/N _  .  -  _ 

HKY-KYl/N 

C _ PARAMETERS  FOR  FET2B  INVERSE  FOURIER  TRANSFORM  SUBROUTINE . 

NRCOEF-N 

_ NCCOEFrN _  _ _ 

LDOOEF-N 

LDA-N 
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C  INITIALIZATION  ROUTINES  FOR  THE  FFT 

_  ....  CALL  FFTCI{N,WFF1)  .  ..  _ _ 

CALL  FFTCI (N,WFF2) 

C„X  IS  THE  WAVENUMBER  IN  DIE  FLUID  ... _ _ 

KH3MEGA/C 

C.  J£C  IS  TOE  COMPLEX  WAVENUMBER  _  _  _ _ 

KC- ( K ) *CMPLX (1 . 0 , < OMEG A*MU* ( ( 4 . 0/3 . 0 ) ♦ ( ( GAMMA-1 . 0 )  /PRA) 

_ i)/(2. 0*RO*C**2. 0) )  )  _  .... _ _ _ 

C  KP  IS  THE  PLATE  WAVENUMBER 

KP“ ( ROB* AQ1B*0M£GA**2 . 0/DB )  *  * . 25  ....  _ 

C  KPF  IS  TOE  COMPLEX  PLATE  WAVENUMBER* *4 . 0 

_  KPF“(ROB*AQlB*CM£Ga**2. 0/DB)*CMPLX(l. O.ETAB)  _ 

DO  100  I-1,N 
DO  100  J"1,N 
BT{ I , J)“0 . 0 

AT(I.J)-0.0  _ 

100  CONTINUE 

C  QMEGAMNSUB  CALCULATES  TOE  UPPER  AND  LOWER  MODES  NUP,  NL0_  ..  . 
CALL  OMEGAMNSUB ( NLO , NUP , OMEGA , ROA , ACHA , DA , LX , LX ) 

C_  TOE  PLATE  MODES  ARE  TAKEN  FROM  EM  AND  EN  -  NLO  TO  NUP  . 

DO  200  EM-NLO , NUP 

_  DO  200  EN-NLO.NUP  _ 

C  OMEGANMSO  IS  TOE  RESONANT  OMEGA* *2 

_  .  QMEGAMNSQ-DSTAR* ( ( ( (2 . 0»EM+T. 0) *PI/LX) *?2 . Q+{ (2»0?EN+1.0). 

**PI/LY)**2.0)**2.0)/(ROA*ACHA) 

_  KXM-( (2.0*EM+1.0)*PI/LX)  _ 

KYN-((2.0*EN+1.0)*PI/LY) 

C.  COEFICIENT  CALCULATES  TOE  INPUT  ARRAY  TO  TOE  FFT  (COEF) _ 

CALL  COEFICIENT  (N,HKX,HKY,KXL,KY1 , PI , KXM,KYN, 

_ tK,KC,KPF(Z,KX,KY,LX,LY,EM,EN,COEF)  _ 

C  FFT2B  IS  AN  IMSL  SUBROUTINE  MUCH  CALCULATES  TOE  2-D  FFT  OF 

C _ A  SET  OF  FOURIER  COEFFICIENTS  _  _  ... _ 

C  COEF ( I / J)  IS  TOE  INPUT  ARRAY 

C_  A(I,J)  IS  TOE  OUTPUT  ARRAY  _ 

CALL  F2T2B  ( NRCXKF , NCCOEF , A , LDA , COEF , LDCOEF , WFF1 , WFF2 , 

_  SWC.CPY)  ...  ..  .  ...... _ 

DO  400  I-1,N 

_ _  DO  400  J-1,N  ....  _ 

B(  I , J)“A( I ,  J) 

_  B1(I,J)-(1.0/(2.0*PI)**2.0)*HKX*HKY*( (-1)**I)*((-1)«J>. 

4*B(I,J)*(COS((2.0»EM+1.0)*PI«XO/LX)*COS((2.0*EN+1.0) 

_  i*PI*YO/LY)/(QMEGAMNSO“QM£GA*«2. 0) )  _  _  _ 

BT(I,J)-BT(I,J)+B1(I,J) 

400  CONTINUE  _ 

200  CONTINUE 

DO  500  I-l.N  a  _ 

DO  500  J-l.N/^  ° 

_  WB(I,J)-((-\.0)*F*RO*QMEGA**2.0*BT(I,J))/(MP*DB) _ 

WBREAL(I.J)-WB<I,J) 

C.  WBREAL(I.J)  IS  TOE  REAL  PART  OF  TOE  PLATE-B  RESPONSE  _ 

C  KB( I , J )  IS  TOC  COMPLEX  PLATE-B  RESPONSE 

500  CONTINUE  _ 

C  TOE  FOLLOWING  TWO  DO-LOOPS  CONVERT  THE  OUTPUT  PLATE  RESPONSE 
C.  IN  TOE  X-Y  PLANE  SO  THAT  TOE  RESPONSE  SHOWS  TOE  CENTER  OF  THE 
C  PLATE  IN  THE  CENTER  OF  TOE  ARRAY 

C  ~  12  112  " 

C.  3  4  3  4  _ 


DO  600  I"1 ,N 
DO  600  J-l.N/2 
TEMP-WBREAL(I.J) 

WBREAL ( I , J ) -WDREAL( I , N/2+ J ) 


WBREAL  ( I ,  N/2+ J  )  -TEMP 

600 _ CONTINUE _  _ _ _ _ 

DO  700  J-l ,  N 

_ DO  700  I-l,N/2 _ _  _ 

TEMP-WBREAL(I,J) 

_ WBREAL(I,J)-WBR£AL(N/2+I,J)  .  ..  . . .  . 

WBR£AL(N/2+I,J)-TEMP 

JOQ _ CONTINUE _ _ _ _ 

AMPB-0.0 

_ COUNT- 0.Q _  _  _  _  .  _  _ 

C  CALCULATE  THE  AMPLITUDE  OF  VIBRATION  IN  TOE  PORTION  OF  TOE  PLATE 

C -SHADOWED. BY  TOE  UPPER  PLATE  ... _ _  _  _ _ 

DO  800  I-(N/2)-<KXl*LX/(4.0»PI)),(N/2)+(KXl»LX/(4.0«PI)) 

- DO  800  J=(N/2).-<KYl«LY/{4.0*PI)),  (N/2)t(KYl*LY/(4.0»EI)> _ 

AMPB- AMPB+ WBREAL ( I , J ) »  »  2 . 0 

_ COUNT-COUNT+l.O- . .  ...  _ 

800  CONTINUE 

C  AMPB  IS  TOE  SUM  OF  TOE  SQUARES  OF  TOE  PLATE  RESPONSE  OF  TOE  AREA 
C  OF  PLATE  B  SHADOWED  BY  PLATE  A 

C _ AMETB .  IS  AMPB  DIVIDED  BY  TOE  NUMBER  OF  POINTS  _ 

AMPTB-AMPB/COUNT 

C_ CALCULATE  THE  VIBRATION  RESPONSE  OF. ELATE  A _  _ 

DO  250  EM-NLO , NUP 

_ DQ-25Q  EN=NLO.NUP. _  ...  .  .  _ 

OMEGAMNSQ-DSTAR*(( ( (2 . 0«EM+1 . 0)«PI/LX) »«2 . 0+( (2.0«EN+1.0) 

_ 4*EI/LY)*?2.Q)**2.Q)/(ROA«ACHA)  _.  .  ..  ...  .  _ 

DO  450  I-1,N 

_ DQ  450. J-l, N _ _  _  _ 

X-(-LX/2.0)+(LX«(I-l))/N 

_ Y-(=LY/2.Q)+(LY«(J-1))/N _ _ _ 

A(I,J)-(COS((2.0«EM+1.0)*PI*XO/LX))«(COS((2.0*EN+1.0) 

- fc*P2*YO/LY) }1(CQS( (2.0*EM+1.0)*PI*X/LX) )!(COS( (2.0*EN+l.fl) _ 

4*PI«Y/LY)  )/(OMEGAMNSO-OMEGA*»2 . 0) 

- AI(I,J)£AT(I,J)+A(I,J)  .  - -  - 

[  450  CONTINUE 

_25Q _ CONTINUE . . .  . 

AMPA-0.0 

_ CQUNT-0.0 _ _ 

DO  550  I-1,N 

_ DQ_55Q  Jsl,N _  _ 

WA(I,J)-(4.0»F/MP)*AT(I,J> 

C  W AREAL  IS  TOE  REALPARTOF  TOE  PLATE  RESPONSE  OF  PLATE  A 

C _ HA— I S-XHE.  COMPLEX  PLATE  RESPONSE  OF  PLATE  A  _  _ 

AMPA-AMPA+WAREAL< I ,  J) ««2 . 0 

_ COUNT-COUNE±l.Q _  _ 

550  CONTINUE 

C  AMPA  IS  TOE  SUM  OF  TOE  SQUARES  OF  THeTpLATE  RESPONSE  OVER  PLATE  A 

C _ AMETA  IS -AMPA  .DIVIDED  BY  TOE  NUMBER  QE  POINTS  .  _ 

WRITE  (26,321)  AMPTA, AMPTB 

C _ OUTPUT  AMPTA  AND  AMPTB _ _ _ _ 

321  FORMAT  <25X,’  AMPTA-  ' ,E12.6,'  AMPTB-  ',£12. 6) 

C_  THE  FOLLOWING  IF  STATEMENT  TESTS  IF  AMPTB  IS  TOO  SMALL  _ 

IF(AMPTB.LE.  (10.0««(-90.0)))GO  TO  650 

_ TAU-AMPTB/AMPTA _  _  _ 

C  CALCULATE  TRANSMISSION  LOSS 

I _ XL(CMEGAI)slQ.Q»(ALOG10(1.0/TAU)) _ _ _ 

Xl(OM£GAI )- FLOAT (OMEGA I } 

C _ OUTPUT  TRANSMISSION  LOSS  __  .  .  . .  . 

WRITE(26,311)OMEGA/(2.0*PI) ,TL(OM£GAI) 

JJJ _ FORMAT (F12. 2, 10X.F6 .2)  .  .  _  _ _ 

GOTO  1000 

.650 _ TL(QMEGAI)?101.0.. .  ...  _  ... 

WRITE(26,311)OMEGA/(2.0*PI) ,TL(OM£GAI) 

1000  CONTINUE 
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SUBROUTINE  COEFICIENT  (N,HKX,HKY,KXl,KYl,PI ,KXM,KYN, 

_ 4K.KC.KPF,Z.KX,KY,LX,LY,£M,£N,C0EF)  . . 

C  COEFICIEf/T  CALCULATES  THE  FOURIER  COEFICIENTS  FOR  USE  IN  THE 

C_  FFT2B  SUBROUTINE-  ..  -  _  .. 

REAL  KX ( N ) , XY  <  N ) ,  J  MN , HKX , HKY , P I , KXI , KYI , 

_  lKXM,KYN,Z,LX,LY,K  ..  _ 

INTEGER  EM,EN,N 

_ COMPLEX  TEMP,KC,KPF,COEF(N,N)  .... _ 

DO  10  I-1,N 

_  DO  10  J-1,N  .  _ 

C  THE  FFT  USES  THE  INTERVAL  FROM  -KXl/2.0  TO  KX1/2.0 

C_  AND  -KYI/2.0  TO  KYl/2.0  ..  _ 

KX(I)-((I-1) *HKX)~( KX1/2 . 0) 

_  KY(J)-((J-l)*HKY)-(KYl/2.0) 

C  THE  FOLLOWING  IF  STATEMENTS  DETERMINE  IMN 

IF  ( (KXM**2.0  .EQ.  KX(  I )  »*2 . 0) .  AND.  _  ....... 

& ( KYN**2 . 0  .EQ.  KY(J>**2.0) ) 

_  &lMN-LX*LY/4.0  _ 

IF  {(KXM**2.0  .EQ.  KX(I)«*2.0) .AND. 

_  fc(KYN««2.0  ,.NE.  KY(J)**2.0))  _ 

1 1MN-LX  *  KYN*COS ( KY { J ) *LY/2 . 0 ) * ( -1 . 0 ) *  * EN/ 

_ 6 ( KYN • • 2 . 0 -KY ( J ) *  *  2 . 0 )  _ _ 

IF  ( (KXM*«2.0  .NE.  KX( I )««2 .0) . AND. 

_ i(KYN**2.0  .EO.  KY( J) **2 . 0) )  .  _ 

iIMN-LY*KXM*COS ( KX ( I ) *LX/2 . 0 ) • ( -1 . 0 ) • *EM/ 

_  A(KXM**2.0-KX(I)**2.0)  .  ...  . . . 

IF  <(KXM**2.0  .NE.  KX< I ) *«2 . 0) . AND. 

_ a{Kyn;«2.0-.ne.  ky(J)«*2.0))  .  _ _ 

i IMN-4 . 0«KXM» KYN*COS ( KX ( I ) • LX/2 . 0 ) * 

_ *COS(KY(J)*LY/2.0)*<t1.0)**EM*<-1.0)**EN/ _ 

t( (KXM**2. 0-KX{ I)**2. 0)* (KYN**2. 0- 

_  tKY( J ) *?2 . 0) )  .  _ 

C  THE  FOLLOWING  IF- STATEMENT  TESTS  IF  THE  EXPRESSION  IS 

C_ NEGATIVE  FOR  USE  IN  THE  COMPLEX  EXXPRESSION  FOR  COEF(I.J) _ 

IF  (K**2.0-KX(I)»*2.0-KY(J) 

_ S”2.0)  55.  65.  65  _ 

55  TEMP-CMPLX(((ABS(K*«2.0-KX(I)*«2.0- 

_ *KY(J)«?2.0))?!.5)*<-Z) ,0.0)  _ _ 

GOTO  75 

_65 _ TEMP“CMPLX(0.0,  <  (K**2.Q-KX(I)**2.0-  .... _ 

SKY(J)«*2.0)«*0.5)*Z) 

-25 _ CONTINUE  .....  _ 

C  CALCULATE  COEF 

_  COEF(I,J)-C£XP{ TEMP)* IMN/  .  _  _ 

A((KC**2.0-KX(I)**2.0-KY<J)**2.0)**.5« 

_ &( (KX{ I)»*2. Q+KY{ J)**2. 0) **2 . 0-KPF) }  _  _ 

10  CONTINUE 

_  RETURN  _ 

END 


SUBROUTINE  OMEGAMNSUB ( NLO . NUP , OMEGA , ROA , ACHA.DA.LX.LY) 

C- . QMEGAMNSUB  COMPUTES  THE  UPPER  AND  LOWER  MODES  WilCli  THE  PLATE 
C  RESPONSE  IS  SUMMED  OVER.  NLO  IS  THE  LOWER  LIMIT  AND  NUP  IS  THE 
C_  UPPER  LIMIT. 

INTEGER  NLO, NUP, Nil 

_ REAL  OMEGA, ROA, AaiA.DA,LX.LY,QMEGAMNO.OMEGAMNl,QM£GAMN2, PI  . 

PI-2. 0*ASIN< 1.0) 

C_  INITIALIZE  AND  INCREMENT  NUP  _ _ 

NLO— 1 

158—  NLD-NLO+1  _ _ 

C  CALCULATE  RESONANT  OMEGA  FOR  THE  EN-NLO,  EM-NLO  MODE 

_  QMEGAMN0-(DA/{ROA*AaU))**.5*(C(2.0*NLOU.0)*PI/LX)**2.0 

*♦( (2 . 0*NLO+1 .0) *PI/LY)**2 .0) 

C  IF  THIS  RESONANT  OMEGA  IS  LARGER  THAN  OMEGA  GO  ON  TO  258 


C  IF  NOT  INCREMENT  NLO  AND  REPEAT 

_ IE.  (OMEGAMNQ.gr. OMEGA)  GOTO  258 _ 

GOTO  158 

C..  INITIALIZE  ANIL  INCREMENT  NUP _ _ 

258  NUP— 1 

358 _ NUPrNUP+1 _ _ _ _ 

C  CALCULATE  THE  RESONANT  OMEGAS  FOR  EN-NUP,  EM-0  AND  EN-O,  EM-NUP 

C  MODES  AND-TEST  IF  LARGER  niAN  OMEGA _  _  _ 

0MEGAMN1-(DA/(R0A«ACHA)  )••  .  5«(  ( (2.0«NUP+1. 0)«PI/LX)»«2. 0 

_ A*<PI/LY);«2.0) _ _  ..  _  _ _  _ 

QMEGAMN2- ( DA/( ROA» ACHA ) ) • • . 5 • { ( P I/LX ) • *2 . 0 

_ a*{<2.o;nue+i.0)*pi/ly>;*2.0)  .  ....  ... _  _ _  ..  _ 

C  IF  LARGER  WAN  OMEGA  GO  ON  TO  458  IF  NOT  INCREMENT  NUP  ,  REPEAT 

- IE(  (QMEGAMNl.GT. OMEGA)  .OR.  (QMEGAMN2.GT. OMEGA) )  GOTO  458 _ 

GOTO  358 

458.  .  CONTINUE . . .  .  .  . 

C  ADD  THREE  MODES  TO  UPPER  LIMIT 

_ NUP-NUP+3 _ _ _ _ _ _ 

Nil -NLO 

C — SUBTRACT  THREE  MODES  FROM  THE. LOWER  LIMIT, .  ZERO  .IS  THE  LOHEST_NLQ._ 
NLO-NLO-3 

_ I£(Nll.£Q.O)  NLO-Q _ _ _ _ _ 

IF(Nll.EO.l)  NLO-O 

_ IF(NU,EQ.2)  NLQsO _ _ _ _ 

WRITE  (6,558)  NLO, NUP 

—558 _ FORMAT.  (215) _ 

RETURN 

_ END _ _ _ 


Sample  Input 


1.000  1.000 
7860. 0007860. 000 
.0000184  343.000 
0.2900  10.0 


0.000  0.0000. 0015870. 001S87 

1.204  0.10000  0.10000  1.4000 

200. 0E9  200. 0E9  402.1  402.1 

10000.0  1  0.100 


.706 

21020.0 


